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Introduction 


Breast  cancer  spread  into  the  lymph  nodes  (LN)  is  the  single  most  deadly 
prognostic  factor  identified  to  date.  Understanding  factors  that  regulate  cancer¬ 
cell  colonization  of  the  lymphatics  is  essential  to  prevent  malignant  invasion. 

Most  breast  cancers  are  estrogen  (ER)  and/or  progesterone  receptor  (PR) 
positive,  and  the  majority  of  breast  cancers  that  spread  to  the  LN  maintain  ER 
and  PR  expression.  However,  due  to  lack  of  experimental  models  of  ER+  breast 
cancer  metastasis,  little  is  known  about  the  roles  of  hormones  in  breast  cancer 
spread  to,  and  growth  in,  LNs. 

Hypotheses:  I  postulated  that  breast  cancer  metastasis  to  LNs  is  dependent  on 
tumor  growth  and  tumor  size.  Second,  I  predicted  that  tumor  spread  to  LNs  is 
dependent  on  tumor  associated  lymphangiogenesis.  Finally,  I  hypothesized  that 
ER  expression  and  transcriptional  function  is  maintained  in  tumors  and  LN 
metastases. 

Project  goals:  I  planned  to  develop  and  then  define  models  of  ER+  breast 
cancer  metastasis,  in  preparation  for  studies  of  hormone  action  that  addressed 
hormonal  regulation  of  key  genes  involved  in  lymphatic  development,  LN  spread, 
and  growth  of  cancer  cells  within  the  LNs. 

Approved  STATEMENT  OF  WORK: 

Task  1.  Develop  and  characterize  all  cell  lines  used  in  the  study  (Months  1- 

6). 

A.  Make  ZsGreen  and  DsRed-Express  polyclonal  stable  cell  lines. 

B.  Quantitation  of  ER  and  PR  protein  expression  in  each  lines. 

Task  2.  Assess  the  impact  of  systemic  estrogen  and  progesterone  on  ER+/- 
PR+/-  breast  cancer  lymphatic  spread  (Months  6-36). 

A.  Grow  and  determine  the  rate  of  lymphatic  invasion  using 
homogeneous  or  dual  colored  mixed  ER+/-,  PR+/-  tumors  with  no 
systemic  hormones,  E,  P,  E+P. 

B.  Perform  immunohistochemistry  on  tumors,  lymphatics,  and  LNs 
to  understand  effects  of  E  and  P  on  cancer  cells,  lymphatics,  and 
LNs. 

Task  3.  Determine  hormonal  influence  to  tumor  lymphangiogenesis 
(Months  8-36). 

A.  Determine  hormone  influence  on  gene  expression  of 
lymphangiogenic  factors. 

B.  Immunohistochemically  quantitate  lymph  vessels  in  each  tumor 
line  and  determine  the  influence  of  hormones  on  lymphatic  vessel 
development. 


4 


Body 


Task  1.  The  majority  of  the  work  for  this  task  is  summarized  in  my  2nd  and  3rd 
first-author  publications  enclosed  within  the  appendix  (1-2)  and  within  the  year  2 
annual  report.  In  these  papers  and  the  previous  report,  I  described  the  studies 
that  complete  and  expand  on  task  1  of  the  approved  statement  of  work.  This 
includes;  (a)  using  fluorescently  labeled  ER+  breast  cancer  cells  (MCF-7  and 
T47D)  to  immunohistochemically  quantify  the  amounts  of  ER  and  PR  expression 
in  the  breast  cancer  cells  grown  in  vivo  (Harrell  et  al,  CanRes  2007  Figure  2). 

This  work  expands  on  my  2006  first  author  publication  (Harrell  et  al,  CanRes 
2006,  figure  6).  By  using  estradiol  withdrawal  post  LN  metastasis  I  was  able  to 
define  the  genes  that  are  estrogen  regulated  in  primary  breast  tumors  compared 
to  their  LN  metastases  (Harrell  et  al,  CanRes  2007  Figure  3.  We  confirmed  that 
ER  were  highly  expressed  in  mammary  gland  tumors  of  all  estradiol-free  mice, 
and  they  behaved  normally,  with  ligand  dependent  downregulation  upon  estradiol 
treatment  coupled  to  PR  upregulation.  However,  these  expression  patterns  were 
altered  in  metastases  within  the  lymphatic  vessels  or  in  LNs,  with  poor  ER 
downregulation  and  decreased  PR  induction.  Thus,  the  lymphatic  environment 
makes  ER+  cancer  cells  estrogen  insensitive.  To  study  this  further,  I  gene 
expression  profiled  estrogen-treated  and  estrogen-withdrawn  tumors  and 
matched  LN  metastases.  For  my  first  study  (1)  the  microarray  analysis  was 
performed  on  laser-captured  MCF-7  cells  from  tumors  and  LN  metastases, 
instead  of  on  whole-tissue.  Interestingly,  the  array  data  confirmed  the 
immunohistochemistry  in  that  fewer  genes  were  E  regulated  in  LNs  as  compared 
to  primary  tumors  (Harrell  et  al,  CanRes  2007  Figure  6).  Furthermore,  some 
genes  that  were  regulated  in  one  direction  by  estradiol  in  the  primary  tumor  were 
regulated  in  the  opposite  direction  by  estradiol  in  the  LN  metastasis  (Harrell  et  al, 
CanRes  2007  Figure  6.  These  findings  are  very  important  as  they  provide  further 
evidence  that  the  LN  microenvironment  promotes  estrogen  resistance.  My 
second  study  from  this  2nd  year  of  funding  contrasted  how  whole-tissue  versus 
laser-captured  material  alters  LN  metastasis  gene  signatures  (2).  Interestingly,  I 
found  that  mouse  RNA,  which  contaminates  whole-organ  microarray  analyses, 
impacts  the  final  data.  The  mouse  hybridizes  to  and/or  dilutes  out  the  human 
RNA,  on  the  human  specific  chips  (Harrell  et  al,  ClinExpMet  Figure  2).  I  found 
that  whole-tissue  and  laser-captured  isolated  RNA  yield  totally  distinct 
“metastasis  signatures”  (Harrell  et  al,  ClinExpMet  Figure  3)  and  I  speculate  that 
the  laser-capture-derived  data  are  more  accurate.  These  studies  complete  task 
1. 

Task  2.  The  majority  of  the  work  to  complete  task  2  was  completed  during  the 
final  year  of  the  funding  period.  I  performed  studies  that  investigated  how 
progesterone  (or  medroxy-progesterone  acetate  (MPA))  alone  or  in  combination 
with  estradiol  affect  ER+  tumor  growth  and  metastasis.  The  findings  of  reduced 
PR  expression  (from  task  1  above)  in  metastases  suggested  that  treating  mice 
with  estradiol  +  progesterone  may  not  affect  cancer  cells  any  differently  once 
they  become  metastatic — as  the  amount  of  PR  that  can  bind  progesterone  is 
diminished.  Though  progestin  treatment  led  to  a  decrease  in  the  percent  of 
tumors  that  generated  LN  metastases,  it  is  obvious  that  in  this  model  system,  if  a 
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ER+  tumor  grows  to  a  certain  size  it  will  metastasize  to  the  local  or  distal  LNs 
(Figure  1).  Further  analyses  of  how  progestin  dose  and  duration  of  treatment 
effect  metastases  is  warranted.  These  studies  complete  task  2. 
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Figure  1:  MCF7+ZsGreen  tumors  were  grown  with  either  E  (n=15),  E+MPA  (n=25). 
E+P  (n=14),  placebo  (n=6),  MPA  (n=7),  or  P  (n=3)  for  12  weeks  in  ovariectomized 
nude  mice.  Tumor  area  and  presence  of  metastases  were  assessed  by  weekly 
whole-body  imaging  and  by  intravital  imaging  at  necropsy.  The  percentage  of  mice 
with  LN  metastasis  is  indicated  at  each  time  point. 


Task  3.  The  third  task  in  my  approved  statement  of  work  is  to  determine 
hormonal  influence  on  tumor  lymphangiogenesis.  In  my  year-one  first-author 
publication  I  show  that  estradiol  dependent  tumor  growth  correlates  with 
intratumoral  lymphangiogenesis  (3);  intratumoral  lymphatic  vessels  were  only 
identified  in  estradiol  treated  MCF-7  tumors.  Peritumoral  LVs  were  found 
surrounding  nearly  all  tumors,  regardless  of  hormone  treatment,  and  likely 
derived  from  the  normal  lymphatic  vasculature  present  in  the  mammary  gland. 

To  further  investigate  hormonal  roles  in  lymphatic  vessel  development,  during 
year  two  of  my  funding,  I  set  up  collaborations  with  Dr.  Kenneth  Korach  at  the 
National  Institutes  of  Environmental  Health  Sciences  and  Dr.  Darryl  Russell  at 
the  University  of  Adelaide  (South  Australia).  Through  RT-PCR  analyses  we  found 
that  ER  alpha  is  required  in  the  ovary  for  proper  lymphatic  remodeling  that  occurs 
after  ovulation,  whereas  ER  beta  has  no  role  in  this  process.  This  is  the  first 
demonstration  that  estrogen  controls  lymphangiogenesis.  This  project  is  being 
expanded  upon  by  my  collaborators  and  we  expect  a  publication  within  the  next 
year.  These  studies  complete  task  3. 


Ongoing  studies.  The  approved  statement  of  work  is  complete. 
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Key  Research  Accomplishments 


•  Developed  a  model  of  fluorescent  estrogen-dependent  LN  metastasis 
(Harrell  et  al.,  CanRes  2006,  Figures  1-3). 

•  Identified  that  tumor  growth  correlates  with  LN  spread  (Harrell  et  al., 
CanRes  2006,  Figure  3). 

•  Found  that  estrogen-dependent  tumor  growth  stimulates  intratumoral 
lymphangiogenesis  (Harrell  et  al.,  CanRes  2006,  Figure  4). 

•  Found  that  the  LN  microenvironment  increases  the  proliferation  rate  of 
cancer  cells  (Harrell  et  al.,  CanRes  2006,  Figure  5). 

•  Found  that  the  LN  microenvironment  increases  ER  expression  and 
decreases  ER  function  (Harrell  et  al.,  CanRes  2006,  Figure  6). 

•  Identified  that  cells  that  express  CD44  are  selected  for  lymphatic 
metastasis  (Harrell  et  al.,  CanRes  2006,  Figure  7). 

•  Identified  that  globally  there  is  less  estrogen  regulation  of  genes  in  LNs  as 
compared  to  primary  tumors  (Harrell  et  al.,  CanRes  2006,  Figure  3). 

•  Identified  specific  genes  that  are  regulated  by  the  tumor 
microenvironment,  by  hormones,  or  by  both  (Harrell  et  al.,  CanRes  2006, 
Figures  3-6). 

•  Determined  that  both  medroxyprogesterone  acetate  (MPA)  and 
progesterone  have  no  influence  on  MCF7  breast  tumor  growth  or 
metastasis  (Figure  1) 


Reportable  Outcomes 

1-  Manuscripts  See  NIH  Biosketch  in  appendices. 

2- Degree  Completion  &  Future  Employment 

My  thesis  was  completed  in  the  summer  of  2007.  Upon  graduation  I 
stayed  in  the  Horwitz  lab  to  complete  this  project  and  to  set-up  my  postdoctoral 
lab.  I  have  since  accepted  an  offer  to  postdoc  with  Dr.  Charles  M.  Perou  at  the 
University  of  North  Carolina,  and  will  start  in  his  lab  towards  the  end  of  2008.  I 
am  currently  a  finalist  for  the  DOD  Era  of  Hope  Postdoctoral  Award. 

Conclusions 


The  three  years  of  funding  for  the  predoctoral  grant  BC050889  were  highly 
successful.  The  research  conducted  yielded  exciting  results  and  further  tested 
the  first  known  reliable  model  of  ER+  breast  cancer  metastasis  in  vivo.  This 
model  and  methods  have  since  led  to  the  establishment  of  collaborations  with 
numerous  investigators  internationally.  Results  found  that  there  are  differences  in 
estrogen  responsiveness  in  breast  tumors  compared  to  the  matched  LN 
metastases,  suggesting  that  aromatase  inhibitors  and  other  estrogen/ER 
mediated  therapies  may  be  differentially  affecting  cancer  cells — dependent  on 
their  tissue  location.  Also,  I  showed  that  laser-capture  microarray  analysis  is 
more  sensitive  and  accurate  than  whole-tissue  analysis.  These  studies  have 
furthered  my  career  and  yielded  opportunities  I  never  thought  would  be  possible. 
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Abstract 

Breast  cancers  commonly  spread  to  lymph  nodes  (LNs).  If  the 
primary  tumors  are  estrogen  receptor  (ER)  and/or  progester¬ 
one  receptor  (PR)  positive,  then  the  likelihood  that  LN 
metastases  express  receptors  exceeds  80%.  However,  due  to 
lack  of  ER+  models,  little  is  known  about  the  role  of  hormones 
in  breast  cancer  spread  or  the  effects  of  the  LN  microenvi¬ 
ronment  on  hormone  responsiveness.  We  have  developed 
metastasis  models  using  ZsGreen  labeled  MCF-7  and  T47D 
human  breast  cancer  cells.  Tumors  are  tracked  in  living  mice 
by  whole-body  imaging,  and  macrometastases  or  micrometa- 
stases  are  detected  by  intravital  imaging  or  fluorescence 
microscopy.  Tumor  growth  is  estrogen  dependent  and 
required  for  intratumoral  lymphangiogenesis.  Seventy-five 
percent  of  all  tumors  and  >95%  of  larger  tumors  generate 
LN  metastases.  Occasionally  more  distant  metastases  are  also 
observed.  “Triads”  of  primary  tumors,  tumor-filled  draining 
lymphatic  vessels,  and  tumor-filled  LNs  from  the  same  mouse 
show  that  (a)  proliferation,  as  measured  by  5-bromo-2r- 
deoxyuridine  uptake,  is  higher  in  the  LN  than  in  the  primary 
tumor.  ( b )  High  ER  levels  are  extensively  down-regulated  by 
estradiol  in  primary  tumors.  However,  there  is  partial  failure 
of  ER  down-regulation  in  LNs  associated  with  (c)  reduced  PR 
expression.  This  suggests  that  ER  are  dysfunctional  in  the  LN 
microenvironment  and  perhaps  hormone  resistant,  (d)  CD44 
is  sparsely  expressed  in  primary  tumor  cells  but  homoge¬ 
neously  overexpressed  in  cells  transiting  the  lymphatics  and 
populating  LNs.  We  hypothesize  that  CD44  expression  targets 
tumor  cells  for  transport  to,  and  uptake  in,  LNs.  If  so,  the 
CD44  pathway  could  be  targeted  therapeutically  to  slow  or 
prevent  LN  metastases.  (Cancer  Res  2006;  66(18):  9308-15) 

Introduction 

Advanced  breast  cancer  kills  >40,000  American  women  each  year 
and  10  times  that  number  worldwide  (1,  2).  Seventy  to  80%  of 
primary  breast  cancers  are  estrogen  receptor  (ER)  and/or 
progesterone  receptor  (PR)  positive  and  considered  to  be  hormone 
responsive  (3).  Importantly,  if  the  primary  tumors  are  ER+,  >80%  of 
lymph  node  (LN)  metastases  and  65%  to  70%  of  distant  metastases 
retain  their  receptors  (4-6).  Despite  this,  little  is  known  about  the 
role,  if  any,  of  estrogens  or  progestins  in  influencing  the  spread  of 
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tumor  cells  from  the  primary  site,  or  their  deposition  and  growth  at 
metastatic  sites. 

This  lack  of  information  is  due,  in  part,  to  a  lack  of  ER+ 
experimental  metastasis  models.  The  three  ER+  human  breast 
cancer  cell  lines  most  widely  used  for  orthotopic  xenograft  studies 
are  MCF-7,  T47D,  and  ZR75  cells  (7).  They  develop  tumors  in 
mammary  glands  of  nude  mice  in  response  to  estradiol 
supplementation  (8,  9).  However,  most  reports  state  that  such 
tumors  are  poorly  invasive  and  rarely,  if  ever,  metastasize  (9). 
Indeed,  most  metastasis  models  do  not  use  these  cells,  relying 
instead  on  ER—  cell  lines  like  MDA-231,  MDA-435  (10)  and  similar 
cells  that  express  putative  aggressiveness  markers  and  do 
metastasize  in  nude  mice.  Other  studies  analyze  only  late  stages 
of  metastases  by  injecting  ER—  tumor  cells  directly  into  the 
circulation  (11).  Alternative  approaches  for  studying  tumor 
aggressiveness  from  orthotopic  sites  with  ER+  breast  cancer  cells 
involve  their  modification  to  overexpress  oncogenes  like  v-Ha-ras 
(12)  or  HER-2/neu  (13),  growth  factors  like  vascular  endothelial 
growth  factor  (14)  or  fibroblast  growth  factor  (15),  or  transcription 
factors  like  Id-1  (16)  or  Fra-1  (17),  in  an  effort  to  enhance  tumor 
progression  while  suppressing  estrogen-dependent  behavior  (18). 
We  considered,  however,  that  if  the  clinical  data  cited  above  are  a 
guide,  it  should  not  be  necessary  to  suppress  hormone  respon¬ 
siveness  to  achieve  a  metastatic  phenotype. 

Observation  of  metastasis  can  be  difficult,  especially  if  few  cells 
are  involved  or  appropriate  target  organs  are  unexplored.  The 
development  of  fluorescent  proteins  to  genetically  tag  living  cells 
has  greatly  increased  the  likelihood  that  metastases  will  be 
visualized  (19).  Indeed,  jellyfish  green  fluorescent  protein  has  been 
widely  used  for  studies  of  ER—  cells  (20).  Recently,  a  set  of  coral  reef 
fluorescent  proteins  were  identified  that  require  no  cofactors  or 
substrates  (21).  These  fluors  have  been  modified  to  produce 
variants  with  even  brighter  fluorescence  and  enhanced  emission 
characteristics  than  green  fluorescent  protein,  and  codon  usage  has 
been  modified  to  optimize  their  mammalian  expression.  Among 
these  fluors,  ZsGreen  has  unique  excitation  and  emission  patterns 
within  the  visible  spectrum. 

This  report  reevaluates  the  metastatic  potential  of  ER+  MCF-7 
and  T47D  human  breast  cancer  cells  in  orthotopic  mouse  tumors 
using  ZsGreen  as  a  sensitive  new  tracking  method.  Cells  stably 
expressing  this  fluor  and  grown  as  xenografts  in  mammary  glands 
of  nude  mice  were  absolutely  estrogen  dependent.  Tumor 
fluorescence  was  stable  and  superficial  LN  metastases  were 
monitored  in  living  mice.  At  necropsy,  stably  fluorescent  metastatic 
cells  were  reliably  detected  microscopically  in  lymphatics  in  local 
and  distant  LNs  and,  occasionally,  in  distant  organs.  If  the  primary 
tumors  were  ER+  and  PR+,  then  metastases  retained  receptors,  but 
their  expression  levels  and  that  of  the  proliferation  marker 
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5-bromo-2'-deoxyuridine  (BrdUrd)  were  modified  by  the  LN  micro¬ 
environment.  In  addition,  the  hyaluronan  receptors,  CD44,  which 
control  chemotactic  behavior  of  tumor  cells  (22),  were  up-regulated 
in  cells  transiting  the  lymphatics  and  retained  in  LNs.  These  data 
suggest  strategies  for  suppressing  LN  metastases  and  indicate  that 
once  in  the  LNs,  tumor  cells  have  reduced  hormonal  sensitivity. 

Materials  and  Methods 

Cell  lines.  MCF-7  cells  were  originally  purchased  from  American  Type 
Culture  Collection  (Manassas,  VA).  Generation  of  T47D  cells  that  express 
one  PR  isoform  was  previously  described  (23),  and  in  the  present  study  PR-B 
expressing  T47D  cells  were  used  exclusively.  The  ZsGreen-Nl  expression 
vector,  PT-67  packaging  cell  line,  and  pLNCX2  retroviral  vector  were 
purchased  from  Clontech/Becton  Dickinson  (Franklin  Lakes,  NJ).  Fluores¬ 
cent  retroviral  vectors  were  generated  by  cloning  the  DNA  coding  sequence 
of  ZsGreen  into  the  pLNCX2  vector  using  the  restriction  enzymes  Hin  dill 
and  AM  and  T4  ligase.  Plasmids  were  then  transfected  into  PT-67 
packaging  cells  with  a  standard  calcium  phosphate  protocol.  Stable 
retrovirus  producing  cells  were  selected  by  2  weeks  of  treatment  with 
500  pg/mL  G418.  MCF-7  or  T47D  cells  were  plated  at  ~  20%  confluence  and 
incubated  in  filtered,  virus-containing  supernatant.  Cells  were  serially 
transduced  two  to  three  times  for  24  hours  each  round,  at  which  point  they 
exhibited  heterogeneous  green  expression.  Cells  were  then  subjected  to 
G418  selection,  followed  by  aseptic  fluorescence  associated  cell  sorting 
(FACS),  to  isolate  a  homogeneously  bright  green  subpopulation.  These  cells 
were  returned  to  culture  and  remained  bright  green  through  multiple 
passage  generations.  Expression  of  ER  and  PR  in  the  selected  cells  was 
confirmed  by  immunohistochemistry  (see  below). 

Xenograft  tumor  growth  and  metastases.  All  animal  procedures  were 
done  under  a  protocol  approved  by  the  University  of  Colorado  Institutional 
Animal  Care  and  Use  Committee.  Ovariectomized  female  athymic  nu/nu 
mice  were  obtained  from  Harlan  Sprague-Dawley  (Indianapolis,  IN)  at  5  to 
6  weeks  of  age.  To  establish  tumors,  animals  were  anesthetized  with  Avertin 
and  injected  into  the  opening  of  the  lactiferous  duct  of  the  abdominal 
mammary  gland  with  1  million  ZsGreen-expressing  MCF-7  or  T47D  cells  in 
100%  Matrigel  (Becton  Dickinson).  Mice  were  also  implanted  with  silastic 
pellets  containing  cellulose  (10  mg)  or  17  (3-estradiol  (2  mg  +  8  mg  cellulose) 
as  previously  described  (24).  Tumor  size  was  recorded  weekly  with  a  digital 
caliper.  Fluorescent  whole-body  imaging  was  done  weekly  (Illumatool  9900, 
Lightools  Research,  Encinitas,  CA).  Two  hours  before  necropsy  100  mg/kg 
BrdUrd  in  PBS  was  injected  i.p.  At  the  end  of  the  study,  mice  were 
euthanized  by  C02  asphyxiation  and  fluorescent  intravital  optical  imaging 
was  done  by  coupling  the  Illumatool  with  an  Olympus  (Melville,  NY)  SZ-61 
dissecting  microscope  and  Olympus  C-5050  digital  camera.  Organs  were 
removed  for  histology,  immunohistochemistry,  and  fluorescence  microscopy. 

Histology,  immunohistochemistry,  and  fluorescence  microscopy. 
Organs  were  fixed  in  4%  paraformaldehyde  overnight,  paraffin  embedded, 
and  cut  into  4-  to  5-pm  sections.  After  high-temperature  antigen  retrieval  in 
citrate  buffer,  sections  were  stained  with  H&E.  Immunohistochemistry  was 
done  on  parallel  sections  for  1  hour  at  room  temperature,  with  primary 
antibodies  directed  against  ER  (6F11;  Vector  Labs,  Burlingame,  CA),  PR 
(1294;  Dako,  Carpinteria,  CA),  BrdUrd  (Becton  Dickinson),  CK18  (Calbio- 
chem,  La  Jolla,  CA),  CK8/CK18  (5D3;  Novocastra,  Newcastle-upon-Tyne, 
United  Kingdom),  LYVE-1  (Upstate,  Lake  Placid,  NY),  Prox-1  (Abeam, 
Cambridge,  MA),  or  CD44  (Ab4;  NeoMarkers/Labvision,  Fremont,  CA). 
Bound  primary  antibodies  were  detected  with  horseradish  peroxidase 
(HRP)-conjugated  goat  anti-mouse  and/or  goat  anti-rabbit  secondary 
antibodies  (Envision  HRP;  Dako).  Sections  were  counterstained  with 
hematoxylin  and  mounted.  For  immunofluorescence  microscopy,  goat 
anti-rabbit  Alexa  Fluor  555  (red)  and  goat  anti-mouse  Alexa  Fluor  488 
(green)  secondary  antibodies  were  used  (Alexa  Fluors;  Invitrogen,  Carlsbad, 
CA).  For  fluorescence  microscopy,  tissues  were  processed  and  sectioned  as 
above,  paraffin  was  removed  with  xylene,  and  then  mounted  with 
fluoromount  G  (EMS,  Hatfield,  PA). 

Statistical  analyses.  Data  from  immunohistochemistry  experiments  were 
collected  by  an  investigator  who  photographed  tissue  sections  and  a  blinded 


investigator  that  counted  cells  from  random  fields  of  stained  sections  using 
Image  Pro  software  version  4.5.0.29  (Media  Cybernetics,  Inc.,  Silver  Spring, 
MD).  Each  field  was  quantified  either  as  a  percentage  of  positively  stained 
cells  or  an  Allred  score  (ref.  25;  percent  of  labeled  cells  added  to  the  intensity 
of  immunolabeling).  Fields  from  the  same  tissue  section  were  averaged  and 
SPSS  software  version  12.0.2  (SPSS,  Inc.,  Chicago,  IL)  was  used  to  analyze  the 
data.  A  Mann-Whitney  U  test  was  used  to  compare  BrdUrd-positive  HRP- 
stained  cells  in  independent,  nonparametric  samples  of  estradiol-  and 
placebo-treated  mouse  tumors.  Paired  t  test  was  used  to  compare  fluorescent- 
stained  cells  positive  for  CK18/BrdUrd,  ER,  or  PR  of  tumor,  LN,  and  associated 
lymphatic  embolus  triad  from  one  mouse  to  those  of  other  triads  from  other 
mice.  A  marginal  homogeneity  test  was  used  to  compare  ER+  and  PR+  HRP 
expression  on  Allred  scores  in  nonparametric  mouse  triad  samples.  For  all 
tests,  P  <  0.05  was  considered  significant. 

Results 

An  in  vivo  model  of  estrogen-dependent  human  breast 
cancer  LN  metastasis.  The  ER+  human  breast  cancer  cell  lines 
MCF-7  and  T47D  are  commonly  used  to  study  the  biology  of  estro¬ 
gen  action  in  vitro.  These  cells  also  grow  as  estrogen-dependent 
xenografts  in  nude  mice  (8,  24,  26),  but  most  published  reports  state 
that  orthotopic  tumors  fail  to  metastasize.  Because,  clinically,  ER+ 
tumors  metastasize  frequently  (4-6),  we  decided  to  reevaluate  the 
metastatic  potential  of  ER+  human  breast  cancer  cells  in  xenografts 
using  sensitive  new  tracking  methods.  To  this  end,  MCF-7  and  T47D 
cells  were  retro virally  transduced  with  ZsGreen,  G418  selected,  and 
the  brightest  green  subpopulation  was  isolated  by  FACS  sorting. 
Confocal  imaging  showed  that  ZsGreen  was  highly  expressed  in 
both  cell  types,  in  the  cytoplasm  and  nuclei  (not  shown).  Intense 
fluorescence  is  maintained  through  multiple  passages,  the  fluor  is 
nontoxic,  and  photobleaching  studies  show  that  ZsGreen,  a  coral 
reef  protein,  is  much  brighter  and  more  photostable  than  green 
fluorescent  protein,  a  jellyfish  protein  (not  shown). 

Visualizing  lymphovascular  tumor  cell  spread.  One  million 
fluorescent  tumor  cells  in  Matrigel  were  injected  into  the  lactiferous 
duct  of  the  abdominal  (#4)  mammary  gland  of  5-  to  6-week-old, 
ovariectomized  nu/nu  mice  supplemented  with  an  estradiol¬ 
releasing  silastic  pellet  (24).  This  injection  route  led  to  formation 
of  well-circumscribed  tumor  cell  nodules  distant  from  the  subiliac 
LN  in  the  gland,  as  shown  by  a  whole  mount  3  days  after  injection 
(Fig.  L4,  top).  In  contrast,  whole  mount  of  a  tumor  at  12  weeks  shows 
the  subiliac  LN,  still  physically  distinct  from  the  primary  tumor  but 
already  highly  involved  with  tumor  cells  (Fig.  1A,  bottom ). 

This  protocol  was  used  to  track  the  spread  of  ER+  tumors 
beyond  the  #4  mammary  gland  using  whole-body  imaging  of  live 
mice  done  weekly.  Figure  IB  shows  an  image  of  tumor  tracks 
photographed  through  the  skin  of  an  intact  mouse  at  12  weeks. 
Bright  green  bilateral  tumors  are  present  in  mammary  gland  #4  at 
the  groin,  and  both  subiliac  mammary  gland  LNs  contain  tumor 
cells.  Also  clearly  visible  are  tumor  cell  tracks  in  the  epigastric 
collecting  lymphatic  vessels  (LV)  draining  the  tumor  cranially  on 
the  right,  as  well  as  tumor  cells  in  the  right  axillary  LN.  Similarly, 
three  distinct  green  fluorescent  clusters  of  tumor  cells  are  faintly 
seen  in  the  left  lymphatics  and  axillary  LN. 

At  sacrifice,  intravital  fluorescence  microscopy  was  used  to 
examine  tumor  spread  to  visceral  tissues  as  shown  in  Fig.  1C 
and  D.  An  axillary  LN  micrometastasis  is  seen  in  one  mouse  at 
12  weeks  (Fig.  1  C,  left )  and  a  fully  involved  LN  metastasis  in 
another  mouse  at  12  weeks  (Fig.  1C,  right).  Figure  1C  also  shows 
the  afferent  LV  full  of  tumor  cells,  with  the  parallel  epigastric  vein 
devoid  of  such  cells.  Tumor  cells  also  spread  medially  to  the 
abdomen.  Shown  in  Fig.  ID  is  early  invasion  into  a  periaortic  LN 
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Figure  1.  Whole-body  and  intravital  images  of 
ZsGreen-expressing  ER+  human  breast  cancer 
cells  and  their  LN  metastases.  A,  excised 
abdominal  mammary  gland  from  a  3-day-old 
MCF-7  tumor,  in  relation  to  the  subiliac  LN 
(top,  dashed  circle),  and  from  a  12-week-old 
MCF-7  tumor  with  a  LN  metastasis  at  a  similar  site 
(bottom).  B,  whole-body  image  showing  bilateral 
T47D  tumors  (T),  subiliac  LN  metastases  (*), 
collecting  LV  (arrows),  and  axillary  LN 
macrometastases  (right)  and  micrometastases 
(left;  dashed  circles).  C,  intravital  image  of  MCF-7 
axillary  LN  metastasis  showing  a  micrometastasis 
(left)  and  a  macrometastasis  (right)  with  a 
collecting  LV  draining  into  it.  D,  intravital  image  of 
an  MCF-7  metastasis  to  the  left  para-aortic  LN 
(left)  and  more  progressive  disease  showing 
coalesced  left  and  right  para-aortic  LNs  and  further 
tumor  spread  caudally  to  the  renal  LNs  (right). 


{left)  and  also  much  more  extensive  metastases  to  periaortic  and 
renal  LNs  {right). 

Bilateral  tumor-filled  lymphatic  tracks  exiting  a  primary  tumor 
are  clearly  seen  in  Fig.  2A.  LV  full  of  tumor  cells  extend  both 
cranially  and  medially  from  the  tumor.  It  is  possible  that  the 
tortuously  branching  network  of  lymphatics  exiting  the  tumor 
cranially  represents  lymphangiogenesis.  A  higher-power  view  of  the 


laterally  branching  iliac  LV  shows  that  the  tumor  cells  within  it 
coalesced  into  emboli  rather  than  moving  as  discrete  cells.  In 
response  to  mechanical  pressure  applied  either  to  the  tumor  or  the 
efferent  lymphatic  root  during  dissections,  fluorescent  emboli  in 
collecting  lymphatics  can  often  be  seen  moving  into  their  draining 
LN  (not  shown).  This  has  implications  for  the  order  in  which  breast 
surgical  procedures  and  sentinel  node  biopsies  are  done.  In  general, 


Figure  2.  Identification  of  cancer  cells  in  LVs. 

A,  top,  intravital  image  of  a  T47D  tumor,  showing 
its  metastasis  to  the  subiliac  LN,  and  cells  moving 
through  the  epigastric  and  iliac  LVs.  Bottom,  higher 
magnification  of  the  iliac  LV  showing  clusters  of 
cells  moving  inwards  toward  the  para-aortic  LNs. 

B,  the  isolated  epigastric  LV  from  (A)  showing 
cancer  cells  moving  through  a  collecting  valve  (*)  in 
route  to  the  axillary  LN.  C,  left,  an  isolated 
collecting  LV  full  of  MCF-7  cancer  cells  (*) 
paralleling  a  milk  duct.  Right,  a  serial  section 
confirms  the  identity  of  the  cancer  cells  and  shows 
maintenance  of  ZsGreen  fluorescence  after  tissue 
processing  into  paraffin  blocks. 
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tumor  cell  spread  to  the  superficial  LNs  is  restricted  to  the  side  of 
the  body  containing  the  primary  tumor,  and  spread  to  deeper  LNs 
becomes  bilateral  with  progressive  disease  (Fig.  ID).  Directional 
valves  {asterisk)  within  the  LV  do  not  inhibit  tumor  cells  moving 
through  LVs,  as  shown  in  Fig.  2 B.  Lastly,  fluorescence  microscopy 
can  be  used  to  analyze  clinicopathologic  features  of  tumors  and 
their  metastases.  Importantly,  ZsGreen  fluorescence  is  maintained 
throughout  the  histologic  fixation  process.  Figure  2C  shows  the 
H&E  stain  of  a  collecting  LV  filled  with  tumor  cells  within  the 
mammary  gland  fat  pad.  A  milk  duct  runs  alongside  it. 
Fluorescence  microscopy  of  a  serial  section  confirms  that  tumor 
cells  are  restricted  to  the  LV. 

Estrogen-dependent  tumor  growth  and  rate  of  LN  meta¬ 
stases.  Figure  3  summarizes  the  primary  tumor  growth  pattern  of 
MCF-7  and  T47D  cells  during  12  weeks  of  observation  and  the 
extent  of  LN  metastasis  as  assessed  by  whole-body  and  intravital 
imaging.  Growth  from  both  cell  types  is  estradiol  dependent,  with 
tumors  nearly  tripling  in  size  in  —  12  weeks  (Fig.  3 A).  Both  MCF-7 
and  T47D  tumors  produce  LN  metastases  at  the  same  rate,  which 
increases  with  time  and  tumor  size.  At  2  weeks,  no  metastases 
occurred.  By  4  weeks,  —  14%  of  tumors  generated  LN  metastases. 
This  increased  to  —75%  at  12  weeks.  Development  of  LN 
metastases  seems  to  be  size  dependent  (Fig.  3 B).  A  separate 
time-course  study  of  tumors  grown  for  1,  6,  or  12  weeks  showed 
that  metastases  to  superficial  subiliac  and  axillary  LN,  or  to  deeper 
visceral  lateral  iliac  and  aortic  lumbar  LN  surrounding  the  aorta 
and  kidneys,  tend  to  occur  at  the  same  rate. 

Estradiol  and  LV  development.  Because  only  estradiol-treated 
tumors  yielded  LN  metastases,  we  next  studied  the  effects  of 
estradiol  on  tumor-associated  LV  development.  For  these  studies, 
tumors  were  grown  with  or  without  estradiol  supplementation  for 
short  term  (3-10  days)  or  long  term  (6-12  weeks)  and  LV  were 
quantified  with  lymphatic-specific  LYVE-1  and  Prox-1  (not  shown) 
antibodies.  Peritumoral  LVs  were  identified  in  both  treatment 
groups  and  likely  represent  lymphatics  found  adjacent  to  veins  and 
arteries  in  the  normal  mouse  mammary  gland  (Fig.  4).  Intratumoral 
LVs  were  not  seen  after  short-term  estradiol  treatment  or  in  the 
estradiol-deprived  dormant  tumors.  Fifty-eight  percent  of  long-term 
estradiol-treated  tumors  contained  intratumoral  LVs  (Fig.  4). 
ZsGreen  fluorescence  was  quenched  during  the  immunolabeling 
process,  allowing  for  double  immunostaining  with  anti-LYVE  and 
anti-BrdUrd  antibodies  (Fig.  4)  to  assess  the  proliferation  state  of 
the  LV.  This  procedure  showed  the  LVs  to  be  BrdUrd  positive, 
suggesting  that  lymphangiogenesis  was  occurring.  However,  double 
labeling  for  LYVE  and  ER  expression  showed  absence  of  ER  in  the 
LVs  (not  shown).  We  conclude  that  LV  proliferation,  although 
restricted  to  the  estradiol-dependent  growing  tumors,  is  either  not 
itself  estradiol  dependent  or  is  occurring  in  response  to  paracrine 
signals  from  surrounding  ER+  tumor  cells. 

The  triad  of  primary  tumor,  lymphatic  emboli,  and  LN 
metastases:  proliferation.  Our  ability  to  dissect  a  primary  tumor, 
its  draining  lymphatic  collecting  duct  filled  with  tumor  cells,  and 
the  downstream  tumor-filled  LN,  from  the  same  mouse,  allows 
assessment  of  the  role  of  the  microenvironment  in  tumor  cell 
proliferation  and  steroid  receptor  expression.  BrdUrd  uptake, 
measured  with  BrdUrd  primary  antibodies  coupled  to  a  HRP  or 
a  green  fluorescent  secondary  antibody,  was  used  to  assess  the 
proliferation  rate  of  cells  in  these  different  compartments.  In 
the  absence  of  estradiol,  tumor  cells  at  the  injection  site  remain 
viable  and  proliferate  slowly  with  a  labeling  index  <7%  (Fig.  5 A), 
but  at  a  rate  insufficient  to  elicit  an  overall  increase  in  tumor 


size  at  12  (Fig.  3)  or  16  to  20  weeks  (not  shown).  These  dormant 
tumors  yield  no  LN  metastases.  In  the  presence  of  estradiol, 
primary  tumors  exhibit  extensive  BrdUrd  uptake  with  a  labeling 
index  of  —  20%. 

To  accurately  quantify  dispersed  tumor  cells  and  eliminate  the 
contribution  of  proliferating  normal  host  cells,  tumor  cells  were 
immunolabeled  with  CK18  and  a  red  secondary  antibody.  In  this 
case,  proliferating  cells  were  quantified  with  anti-BrdUrd  and  a 
green  fluorescent  secondary  antibody,  and  some  sections  were 
also  counterstained  with  blue  4r,6-diamidino-2-phenylindole 
(DAPI).  Figure  5 B  shows  two  “triads”  of  tumor,  LV  and  LN.  Each 
triad  came  from  the  same  mouse;  the  two  sets  were  taken  from 
two  different  mice.  In  triad  1,  although  the  primary  tumor  cells 
have  a  high  labeling  index,  proliferation  in  the  LV  embolus  is  low 
and  mainly  restricted  to  cells  at  the  periphery,  and  in  the  LN, 
tumor  cells  are  present  only  in  the  subcapsular  sinus.  The 
remainder  consists  of  normal  DAPI-stained  LN  cells  (some  of 


Figure  3.  Estrogen-dependent  tumor  growth  is  required  for  LN  metastasis. 

A,  MCF-7  ( triangles )  and  T47D  ( circles )  tumors  were  grown  either  with  ( +E )  or 
without  (-E)  estradiol  for  12  weeks  in  ovariectomized  nude  mice.  Tumor  area 
was  quantified  with  digital  calipers  and  presence  of  metastases  was  assessed 
by  weekly  whole-body  imaging  and  by  intravital  imaging  at  necropsy.  The 
percent  of  mice  with  LN  metastasis  is  indicated  at  each  time  point.  B,  tumor  area 
at  necropsy  as  assessed  by  the  diameter  in  two  dimensions.  Filled  symbols, 
tumors  from  which  metastasis  developed;  open  symbols,  tumors  that  did  not 
generate  metastases. 
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Figure  4.  Assessment  of  tumor-associated  lymphatics  and  lymphangiogenesis. 
Tumors  from  3-  to  10-day-old  or  6-  to  12-week-old  placebo-treated  (-E) 
and  estradiol-treated  ( +E )  ovariectomized  mice  were  paraffin  embedded. 

Top,  sections  (5  |am)  were  probed  for  LVs  with  an  anti-LYVE-1  antibody  and  an 
HRP-labeled  secondary  antibody.  Bottom,  sections  were  probed  by  dual 
fluorescence  immunohistochemistry  with  mouse  anti-Brdllrd  and  rabbit 
anti-LYVE  primary  antibodies,  followed  by  goat  anti-mouse  (green)  and  goat 
anti-rabbit  (red)  secondary  antibodies,  to  identify  lymphatics  in  S  phase.  The 
table  quantifies  the  results. 


which  are  also  undergoing  cell  division).  Triad  2  exhibits  a 
different  pattern,  with  extensive  proliferation  in  the  primary 
tumor  and  throughout  the  LV  emboli,  and  the  LN  is  highly 
involved  with  rapidly  proliferating  tumor  cells. 

Figure  5 C  quantifies  all  proliferation  data.  It  shows,  first,  that  the 
two  methods  of  assessing  BrdUrd  uptake  yielded  similar  results 
(compare  the  two  +E  “tumor”  sets— one  assessed  with  HRP  and  the 
other  with  dual  BrdUrd/ CK18  immunohistochemistry).  As  expected, 
compared  with  dormant  controls,  estradiol-treated  primary  tumors 
exhibit  a  statistically  significant  increase  in  proliferation  rate  from 
~  7  to  ~  20%.  And,  compared  with  the  primary  tumors,  metastatic 
cells  in  the  LN  microenvironment  have  a  statistically  significant 
(P  =  0.003)  further  increase  to  ~  27%.  Proliferation  rates  of  embolic 
cells  in  the  LV  were  variable  and  not  statistically  different  from  the 
other  two  compartments,  although  they  tended  to  resemble  those  in 
the  primary  tumor  more  closely. 

The  triad  of  primary  tumor,  lymphatic  emboli,  and  LN 

metastases:  ER  and  PR.  ER  and  PR  immunohistochemistries  for 
two  sets  of  vehicle  versus  estradiol-treated  triads  using  MCF-7 
derived  tumors  show  the  heterogeneity  observed  (Fig.  6).  In  the 
absence  of  estradiol,  levels  of  ER  in  tumors  are  always  high. 
Estradiol  treatment  leads  to  extensive  tumor  ER  down-regulation. 
ER  tend  to  remain  low  in  the  LV  and  LN  of  Fig.  6  {top)  but  are 
somewhat  restored  in  the  LN  of  the  second  triad  {bottom).  No  PR 
are  expressed  without  estradiol  induction.  Estradiol  treatment 
leads  to  heterogeneous  PR  expression  in  tumors,  which  tends  to 
decrease  in  LV  and  LN  of  the  same  mouse.  Quantitation  of  PR 


data  from  several  mice  is  shown  in  Fig.  6.  There  is  a  clear  trend 
of  decreased  PR  expression  as  cells  spread  into  LV  and  LN 
{P  =  0.052).  A  second  experiment  (not  shown)  using  fluorescent 
immunohistochemistry  for  ER  and  PR  showed  a  similar  trend 
of  reduced  PR  between  five  tumors  and  their  matched  LNs 
{P  =  0.09).  Each  immunohistochemical  procedure  showed  a 
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Figure  5.  Assessment  of  tumor  proliferation  in  MCF-7  tumors,  LV  emboli, 
and  LN  metastases.  A,  ovariectomized  mice  bearing  MCF-7  tumors  were  treated 
for  12  weeks  with  placebo  (-E)  or  estradiol  (+E).  Tumors  were  paraffin 
embedded,  sectioned,  and  stained  with  H&E  or  analyzed  by  BrdUrd-HRP 
immunohistochemistry.  B,  dual  fluorescence  immunohistochemistry  for 
cytokeratin  18  (red)  and  BrdUrd  (green)  of  estradiol-treated  MCF-7  tumors, 

LV  emboli,  and  LN  metastases.  Each  triad  was  taken  from  the  same  mouse. 
Anti-rabbit  CK18  and  anti-mouse  BrdUrd  primary  antibodies  and  goat  anti-rabbit 
(red)  and  goat  anti-mouse  (green)  antibodies  were  used.  A  DAPI  counterstain 
(blue)  was  occasionally  included.  The  centrally  necrotic  regions  of  tumors  were 
excluded  from  analyses.  C,  quantification  of  A  and  B  (n  =  5  for  tumors  and  LNs; 
n  =  3  for  LVs),  showing  the  data  obtained  from  the  HRP  based  protocol  in 
A  (white  columns)  and  the  dual  CK18/BrdUrd  protocol  in  B  (black  columns). 
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Figure  6.  ER  and  PR  expression  in  triads  of  MCF-7  tumors,  LV  emboli,  and  LN  metastases.  Left,  representative  images  of  two  tumors  from  placebo  (-E)  treated 
mice  and  of  two  matched  triads  of  estradiol-treated  ( +E )  mice,  showing  tumors,  their  lymphatic  emboli,  and  LN  metastases.  Paraffin-embedded  sections  were 
probed  with  mouse  anti-ER  or  mouse  anti-PR  primary  antibodies,  coupled  with  goat  anti-mouse  HRP-conjugated  secondary  antibodies.  Right,  quantification  of  PR 
expression  in  placebo-treated  (n  =  3)  and  estradiol-treated  tumors  (n  =  7),  LV  emboli  ( n  =  3),  and  LN  metastases  ( n  =  5),  based  on  the  intensity  of  cellular  staining 
combined  with  the  percent  of  positive  cells. 


nonsignificant  slight  increase  in  ER  expression  in  the  LN 
metastases  compared  with  the  down-regulated  ER  in  the 
estradiol-treated  primary  tumor. 

Selection  of  cells  for  lymphatic  transit:  CD44.  Hyaluronan 
regulates  various  aspects  of  cell  behavior,  particularly  cell 
migration  and  invasiveness.  It  does  so  by  binding  to  specific 
receptors  on  the  cell  surface,  the  best  characterized  of  which  is 
CD44  (22).  Human  breast  cancer  cells  migrate  directionally  in 
response  to  a  hyaluronan  gradient  depending  on  CD44  expression 
levels  (22).  Breast  cancers  secrete  factors  that  up-regulate 
hyaluronan  expression  in  osteoblasts  (27),  which  may  promote 
bone  metastasis.  Figure  7 A  shows  CD44  expression  levels  in  two 
primary  tumors  and  in  their  LV  and/or  LN  metastases.  Set  A  (top) 
clearly  shows  extensive  overexpression  of  CD44  in  tumor  cells 
occupying  the  LN  microenvironment,  compared  with  the  same 
tumor  cells  in  the  mammary  gland,  where  CD44  expression  is 
sparse  and  heterogeneous.  Set  B  includes  a  matched  LV  embolus 
and  again  shows  sparse  and  heterogeneous  CD44  in  the  tumor,  but 
homogenous  and  elevated  CD44  in  both  the  LV  embolic  tumor 
cells,  and  the  LN  metastases. 

Distant  metastases.  Metastases  to  distant  organs  were  occa¬ 
sionally  observed  with  these  models.  Figure  IB  (left)  shows  H&E 
staining  of  tumor  cells  metastatic  to  the  lung.  Their  identity  is 
confirmed  by  intense  ZsGreen  fluorescence  in  a  serial  section 
(Fig.  7 B,  right).  These  lung  metastases  could  have  arisen  by  direct 
peritumoral  venous  invasion  and  spread  (the  most  direct  route  to 


the  lungs  from  distant  body  sites)  or  they  may  reflect  indirect 
spread  from  an  overwhelmed  lymphatic  system  via  lymphatic- 
venous  anastomoses.  Both  mechanisms  may  contribute  to  distant 
organ  metastasis  in  humans  (28).  It  should  be  stressed  that  these 
metastases  were  rare  and  they  occurred  only  in  animals  that  also 
had  extensive  LN  spread. 

Perineurial  invasion.  Perineurial  invasion  was  also  noted  in 
tumors  arising  from  both  cell  lines.  Perineurial  invasion  is  observed 
in  a  variety  of  human  carcinomas  and  may  contribute  to  local- 
regional  tumor  spread,  although  it  has  not  been  shown  to  be  an 
independent  prognostic  indicator  in  human  mammary  carcinomas 
(29).  Although  not  as  frequent  as  lymphatic  invasion,  perineurial 
invasion  was  occasionally  found  at  tumor  margins  and  along  the 
peripheral  nerves  that  parallel  the  iliac  veins.  Figure  7C  shows  an 
H&E  stain  (left)  and  CK18  immunohistochemistry  (right)  of  T47D 
tumor  cells  (arrow)  tracking  along  a  peripheral  nerve  (asterisk). 

Discussion 

LN  metastasis  in  association  with  tumor  size  is  the  single  most 
powerful  indicator  of  poor  prognosis  in  mammary  carcinoma 
(30,  31).  Pathologic  data  sets  of  patients  with  breast  cancer 
indicate  that,  at  diagnosis,  presence  of  tumor-infiltrated  LNs 
is  common,  with  estimates  ranging  from  30%  to  50%  of  cases 
(32-34),  depending  on  tumor  size.  The  dominant  feature  of  the 
present  model  is  rapid  and  reliable  LN  metastasis  of  ER+  cells 
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from  orthotopic  tumor  sites.  Two  different  cell  lines,  MCF-7  and 
T47D,  generated  estrogen-dependent  solid  tumors  and  reliably 
produced  metastases  to  local  and  distant  LNs.  Under  conditions 
of  the  studies,  which  were  terminated  at  12  weeks,  tumor  cells  in 
LV  and  LN  metastases  were  observed  in  75%  of  all  mice  and  in 
essentially  100%  of  mice,  the  primary  tumors  of  which  exceeded 
—  75  mm2  in  area  (Fig.  3).  Clinically,  the  association  between 
tumor  size  and  LN  involvement  is  well  known  (35).  Metastases  to 
other  sites,  like  the  lungs  (Fig.  7B),  were  occasionally  observed,  as 
were  cells  migrating  along  peripheral  nerve  tracks  (Fig.  7C).  The 
ability  of  the  latter  routes  to  reliably  colonize  distant  organs 
undoubtedly  will  require  longer  observation  times.  This  is 
suggested  by  the  studies  of  others.  Half  a  million  ER—  MDA-231 
breast  cancer  cells,  injected  directly  into  the  circulation, 
nevertheless  required  10  to  12  weeks  before  exhibiting  limb  bone 
metastases  in  30%  of  mice  (36).  ER+  ZR75  cells  injected  into  the 
circulation  of  mice  required  3  to  6  months  to  yield  similar  bone 
metastases  (37). 

Use  of  modified  breast  cancer  cells  overexpressing  hormones  (38), 
growth  factors  (14,  39),  and  oncogenes  (12)  shows  that  the 


Figure  7.  Selection  of  CD44+  MCF-7  tumor  cells  for  lymphatic  transit  and 
examples  of  lung  metastases  and  perineurial  invasion.  A,  tumors  (n  =  5)  taken 
from  estradiol-treated  mice,  and  their  matched  LV  (r?  =  3)  and/or  LNs  (n  =  6), 
were  probed  with  mouse  anti-CD44  and  goat  anti-mouse  HRP-conjugated 
secondary  antibodies.  Top,  a  tumor  and  its  matching  LN.  Bottom,  a  triad  of 
tumor,  LV,  and  LN  taken  from  the  same  mouse.  B,  left,  H&E  stain  of  a 
paraffin-embedded  section  of  T47D  cells  from  an  estradiol-treated  mouse  with 
metastasis  to  the  lung.  Right,  a  serial  section  processed  for  fluorescence 
microscopy.  C,  left,  H&E  stain  of  MCF-7  tumor  cells  from  an  estradiol-treated 
mouse  showing  perineurial  invasion.  Right,  CK18  stain  of  a  deeper  section. 

*,  nerve;  arrow,  tumor  cells. 


metastatic  potential,  route  of  dispersal,  time  to  metastasis,  and 
preferred  organ  arrest  site  of  any  cell  line  are  dependent  on  unique 
factors  elaborated  by  that  cell,  and  that  much  remains 
to  be  learned  about  the  role  of  estrogens  and  progestins  on 
these  processes.  For  example,  it  has  been  reported  that  clinically,  ER+ 
tumors  are  more  likely  to  spread  to  bones  whereas  ER—  tumors 
more  commonly  metastasize  to  the  viscera  (40,  41).  ER+  metastasis 
models,  reflecting  the  clinical  situation,  could  help  define  molecular 
mechanisms  involved  in  such  tissue  specificity. 

Expression  profiling  of  matched  primary  and  metastatic  LN 
tumor  pairs  in  clinical  samples  has  failed  to  identify  genes  that 
distinguish  one  from  the  other  (42)  because  metastases  are 
genetically  more  closely  related  to  the  primary  tumor  from  which 
they  arose  than  to  any  other  metastasis.  However,  within  matched 
pairs,  select  genes  are  differentially  expressed  between  the  tumor 
and  LN  metastasis.  Common  among  these  are  extracellular  matrix 
and  cell-matrix  interacting  genes,  like  the  hyaluronan  receptor 
CD44,  which  mediates  rolling,  attachment,  and  migration  of  cells 
on  a  hyaluronan  substratum  (43,  44).  We  show  here  that  CD44  is 
sparsely  and  heterogeneously  expressed  in  primary  tumors  but 
homogeneously  overexpressed  in  LV  tumor  cell  emboli  and  LN 
metastases  (Fig.  7).  Because  directional  migration  of  breast  cancer 
cells  towards  hyaluronan  is  dependent  on  CD44  (22),  we  propose 
that  the  CD44+  subpopulation  of  primary  tumor  cells  is 
preferentially  propelled  into  tumoral  lymphatics,  and  from  there 
to  LNs.  If  so,  therapeutic  targeting  of  the  CD44+  subpopulation  in  a 
primary  tumor  could  prevent  LN  metastases. 

There  is  a  high  concordance  rate  (  —  80-90%)  in  ER  expression 
between  primary  tumors  and  matched  LN  metastases  (5).  However, 
in  the  same  patient,  considerable  disparities  can  exist  in  the 
concentrations  of  ER  at  the  two  sites,  with  ER  levels  often  (6,  45)  but 
not  always  (45)  higher  in  the  LN  than  in  the  local  primary  tumor.  We 
observe  a  pattern  (Fig.  6)  in  which  ER  levels  are  higher  in  the  LN 
than  in  the  primary  tumor  of  the  same  mouse,  an  indicator  of  ER 
down-regulation  failure  in  the  LN  microenvironment.  It  has  been 
reported  that  there  is  an  association  between  estrogen-dependent  ER 
down-regulation  and  improved  transcriptional  efficacy  of  the 
receptors  (46).  Although  it  might  be  counterintuitive  to  think  that 
in  the  presence  of  estrogen,  lower  ER  levels  signal  a  stronger 
hormone  response,  this  seems  to  be  the  case.  Thus,  the  inefficient  ER 
down-regulation  in  the  LN  could  be  associated  with  poorer  hormone 
responsiveness  at  that  site.  In  the  mouse  tumors,  PR  expression 
levels  (Fig.  6),  which  are  markers  of  estrogen  action  and  ER  function, 
support  this  hypothesis.  Assessment  of  PR  levels  in  LNs  and  matched 
primary  tumors  shows  a  trend  towards  a  decrease  in  LN  PR  (Fig.  6). 
This  would  be  predicted  by  inefficient  ER  down-regulation.  Clinically, 
both  up-regulation  and  down-regulation  of  PR  have  been  reported  in 
the  LN  compared  with  primary  tumors  (45),  but  this  has  not  been 
correlated  to  ER  in  the  same  samples. 

We  also  observe  a  statistically  significant  increase  in  the 
proliferation  rate  of  LN  metastases  compared  with  matched 
primary  tumors  (Fig.  5 B  and  C).  Clinically,  increased  mitoses  in 
nodal  metastases  compared  with  primary  tumors  have  been 
reported  (47),  possibly  due  to  release  of  cytokines  and  growth 
factors,  such  as  insulin-like  growth  factor-I  and  epidermal  growth 
factor,  by  the  LN  (48).  The  implications  of  this  growth  promotion 
within  LNs  for  failure  of  hormone  treatments  in  metastatic 
disease  could  be  important.  If,  in  the  LN  microenvironment, 
cytokines  and  growth  factors  are  the  dominant  mitogens,  then  LN 
metastases  would  be  less  sensitive  than  primary  tumors  to 
estradiol-suppressive  therapies.  Indeed,  in  metastatic  disease, 
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tamoxifen,  which  targets  ER  signaling  pathways,  is  rarely 
associated  with  long-term  remission,  exhibiting  decreased  mean 
time  to  disease  progression  and  decreased  duration  of  response 
than  local  disease  (41,  49,  50). 

In  summary,  we  have  developed  models  of  ER+  human  breast 
cancers  that  reliably  spread  to  LNs  from  solid  tumors.  Through  the 
use  of  a  new  generation  fluorescent  protein,  ZsGreen,  we  are  able 
to  identify  cancer  cells  in  tumors,  LV  emboli,  LNs,  and,  occasionally, 
more  distant  sites  for  further  molecular  analysis.  These  studies 
show  that  the  tumor  cells  found  in  LNs  overexpress  markers 
associated  with  tumor  aggressiveness.  At  the  same  time,  the  LN 
microenvironment  may  suppress  estrogen  sensitivity  of  tumor 


cells,  which  would  make  them  more  resistant  than  primary  tumors 
to  ER-targeted  treatments. 
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The  h  igh  propensity  of  breast  cancer  to  metastasize  to  bone 
results  from  contributions  on  both  sides  of  the  tumor-stroma 
interface.  Identifying  the  components  in  bone  (growth  factors, 
etc.)  essential  for  metastasis  is  challenging.  To  investigate 
bone-breast  cancer  interactions,  Moreau  and  colleagues 
created  engineered  bone  by  adding  human  marrow-derived 
mesenchymal  stem  cells  to  silk  fiber  scaffolds.  Later, 
metastases  were  demonstrated  to  form  in  engineered  bone 
from  human  breast  cancer  cells  originating  at  the  orthotopic 
site  in  a  mouse  model.  Since  engineered  bone  can  be 
manipulated  in  controlled  fashion,  it  is  now  possible  to  use 
it  to  identify  bone  components  essential  for  metastasis. 

Pregnane  X  Receptor  Enhances 
Chemotherapy  in  Prostate  Cancer 

Chen  et  al. _ Page  10361 

As  a  sensor  of  xenobiotics,  pregnane  X  receptor  (PXR) 
regulates  the  gene  expression  of  drug  metabolizing  enzymes 
and  efflux  transporters  that  can  have  profound  effects  on  the 
efficacy  of  chemotherapy.  Chen  and  colleagues  found  that 
prostate  cancer  cells  expressed  functional  PXR  which,  when 
activated,  rendered  tumor  cells  with  increased  resistance  to 
Taxol  and  vinblastine.  Ablation  of  PXR  in  prostate  cancer  cells 
enhanced  the  efficacy  of  chemotherapy.  Since  PXR  can  be 
activated  by  a  range  of  chemotherapeutics  and  herbal 
supplements,  potential  activation  of  PXR  should  be  taken  into 
consideration  when  designing  combination  chemotherapy 
for  prostate  cancer. 

Mutagenesis  Screen  Indentifies 
EGFR  Inhibitors 

Yu  et  al.  _ _ Page  10417 

Secondary  EGFR  mutations  commonly 
lead  to  resistance  to  the  EGFR 
inhibitors  erlotinib  and  gefitinib  in 
patients  with  non-small  cell  lung 
cancer.  Irreversible  EGFR  inhibitors, 
such  as  CL-387,785  can  overcome 
resistance  and  are  in  clinical  development.  Yu  and  colleagues 
developed  a  cell-based,  in  vitro  random  mutagenesis  screen  to 
identify  EGFR  mutations  that  confer  resistance  to  CL-387,785 


and  identified  a  number  of  novel  resistance  mutations. 

Both  an  alternative  EGFR  inhibitor  and  a  Cdk4  inhibitor 
could  overcome  such  resistance.  These  results  identify 
novel  mutations  mediating  resistance  to  irreversible  EGFR 
inhibitors  and  reveal  alternative  strategies  to  overcome  or 
prevent  resistance  in  EGFR-mutant  lung  cancers. 

Aromatase  Levels  Predict  Survival  in 
Women  with  Lung  Cancer 

Mah  et  al. _ Page  10484 

Lung  cancer  remains  the  most  lethal  type 
of  human  malignancy,  with  increased 
incidences  occurring  especially  in  women. 
Mah  and  coauthors  describe  novel  findings 
about  the  prognostic  value  of  aromatase 
expression  in  predicting  the  course  of 
non-small  cell  lung  cancer  (NSCLC).  This  group 
demonstrated  the  importance  of  estrogen  and  aromatase 
for  promoting  tumor  progression  in  a  mouse  model.  They 
further  conducted  a  population-based  study  examining 
aromatase  levels  in  human  NSCLC  using  tissue  microarray 
technology.  Lower  levels  of  aromatase  were  a  strong  predictor 
of  survival  in  women  older  than  65  years  with  earlier  stage 
cancer.  The  authors  predict  that  the  use  of  aromatase 
inhibitors  may  benefit  this  patient  population. 

Lymph  Node  Microenvironment  Alters 
Sensitivity  to  Estradiol 

Harrell  et  al  _ page  1 0582 

Estrogen-responsive  breast  tumors  often  become  refractory 
to  antiestrogens  and  recur  in  organs  distant  from  the 
primary  site.  To  determine  if  the  lymph  node  (LN) 
microenvironment  alters  tumor-cell  sensitivity  to  estradiol, 
Harrell  and  colleagues  developed  a  model  to  assess  the 
hormone  responsiveness  of  xenografted  human  breast 
tumors  and  their  matched  LN  metastases.  Laser-capture 
microdissection  plus  gene  expression  profiling  confirmed 
by  immunohistochemistry,  found  that  LN  metastases  were 
estrogen  insensitive  compared  to  the  primary  tumors. 

Indeed,  some  genes  were  regulated  by  estradiol  in  the 
opposite  direction  in  tumors  and  matched  metastases. 

This  study  shows  that  the  LN  microenvironment  influences 
expression  of  therapeutic  targets  that  may  alter  cancer  cell 
survival,  recurrence,  or  response  to  treatments. 
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Abstract 

The  lymphatic  system  is  a  common  avenue  for  the  spread  of 
breast  cancer  cells  and  dissemination  through  it  occurs  at 
least  as  frequently  as  hematogenous  metastasis.  Approximate¬ 
ly  75%  of  primary  breast  cancers  are  estrogen  receptor  (ER) 
positive  and  the  majority  of  these  maintain  receptor  expres¬ 
sion  as  lymph  node  (LN)  metastases.  However,  it  is  unknown  if 
ER  function  is  equivalent  in  cancer  cells  growing  in  the  breast 
and  in  the  LNs.  We  have  developed  a  model  to  assess  estrogen 
responsiveness  in  ER+  breast  tumors  and  LN  metastases. 
Fluorescent  ER+  MCF-7  tumors  were  grown  in  ovariectomized 
nude  mice  supplemented  with  estradiol.  Once  axillary  LN 
metastasis  arose,  estradiol  was  withdrawn  (EWD),  for  1  or 
4  weeks,  or  continued,  to  assess  estradiol  responsiveness. 
On  EWD,  proliferation  rates  fell  similarly  in  tumors  and  LN 
metastases.  However,  estradiol-dependent  ER  down-regulation 
and  progesterone  receptor  induction  were  deficient  in  LN 
metastases,  indicating  that  ER-dependent  transcriptional 
function  was  altered  in  the  LN.  Cancer  cells  from  estradiol- 
treated  and  EWD  primary  tumors  and  matched  LN  metastases 
were  isolated  by  laser  capture  microdissection.  Global  gene 
expression  profiling  identified  transcripts  that  were  regulated 
by  the  tissue  microenvironment,  by  hormones,  or  by  both. 
Interestingly,  numerous  genes  that  were  estradiol  regulated 
in  tumors  lost  estradiol  sensitivity  or  were  regulated  in  the 
opposite  direction  by  estradiol  in  LN  metastases.  We  propose 
that  the  LN  microenvironment  alters  estradiol  signaling  and 
may  contribute  to  local  antiestrogen  resistance.  [Cancer  Res 
2007;67(21):10582-91] 

Introduction 

Interactions  between  luminal  epithelial  cells  and  the  surrounding 
microenvironment  govern  the  overall  physiology  of  the  mammary 
gland  (1),  with  development  of  epithelial  ducts  being  dependent 
on  the  activity  of  both  stromal  and  epithelial  estrogen  receptor  a 
(ERa;  ref.  2).  Integrins,  growth  factors,  and  steroid  hormone  sig¬ 
naling  pathways  all  play  an  important  part  in  maintaining  normal 
glandular  architecture  (3).  Stroma  makes  up  >80%  of  breast  volume 
and  is  composed  of  fat,  interstitial/interlobular  dense  connective 
tissue,  intralobular  loose  connective  tissue,  and  lymphatic/blood 
vessels  (4).  Luminal  epithelial  cells  associate  with  these  stromal 
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elements  as  well  as  with  basement  membrane  and  myoepithelial 
cells.  Disruption  of  this  delicate  balance  results  in  dramatic  changes 
in  both  extracellular  and  intracellular  signaling  (3). 

Seventy  percent  to  80%  of  primary  breast  tumors  express  ERa 
and  initially  respond  to  estradiol  withdrawal  (EWD)  or  antiestrogen 
therapies  (5).  Breast  cancer  ERa  expression  and  function  are  also 
influenced  by  the  cellular  microenvironment.  In  vitro  studies 
with  human  ER+  breast  cancer  cells  show  that  extracellular  matrix 
proteins,  such  as  type  I  collagen  and  laminin,  can  modify  proli¬ 
ferative  responsiveness  to  estrogen  (6).  A  comparison  of  in  vitro 
and  in  vivo  models  found  that  estradiol  regulates  different  genes 
in  human  breast  tumor  xenografts  compared  with  the  identical 
cells  in  culture  (7).  These  studies  suggest  that  a  microenvironment- 
dependent  influence  on  estradiol-dependent  gene  expression  exists 
in  ER+  breast  cancers.  However,  this  issue  has  received  little  atten¬ 
tion,  and  to  the  best  of  our  knowledge,  in  vivo  experiments  that  test 
if  metastatic  microenvironments  alter  estradiol-dependent  gene 
expression  have  not  been  conducted. 

Metastasis  to  the  lymph  nodes  (LN)  is  a  key  prognostic  factor 
that  conveys  advanced  disease  status  with  the  possibility  that 
cancer  cells  have  spread  to  other  more  distant  sites.  At  diagnosis, 
30%  to  50%  of  all  breast  cancers  have  spread  to  the  draining  or 
“sentinel”  LN  (8-10),  and  if  the  primary  tumor  is  ER+,  then  ~  80% 
of  LN  metastases  retain  ER  expression  (11-13).  As  tumor  cells 
disseminate  to  LNs,  they  enter  a  distinct  microenvironment.  Here, 
they  encounter  different  supporting  elements,  such  as  reticular 
fibers,  flowing  lymph  fluid,  lymphocytes,  T  and  B  cells,  increased 
numbers  of  macrophages,  and  epithelial  reticular  cells  that  are  sur¬ 
rounded  by  adipose  tissue  (14).  It  is  unknown  if  the  differences  in 
structure,  function,  and  composition  of  the  mammary  gland  com¬ 
pared  with  the  LN  affect  ER  activities  and  estradiol-dependent  gene 
expression  in  breast  cancer  cells  that  are  located  within  these  two 
distinct  microenvironments. 

The  studies  presented  here  use  a  xenograft  model  of  fluorescent, 
ER+,  estradiol-dependent  human  breast  cancer  LN  metastasis  to 
address  this  question.  In  the  same  mouse,  LN  metastases  are  found 
to  be  transcriptionally  estradiol  insensitive  compared  with  primary 
tumors  from  which  they  originated,  despite  retention  of  ER.  We 
identify  genes  that  may  aid  in  explaining  mechanistically  why 
advanced  ER+  breast  cancers  become  nonresponsive  to  antiestro¬ 
gen  or  EWD  therapies  and  often  recur  as  metastases. 

Materials  and  Methods 

Cell  lines.  MCF-7  human  breast  cancer  cells  were  originally  purchased 
from  the  American  Type  Culture  Collection.  The  generation  of  fluorescent 
MCF-7  cells  has  been  previously  described  (15).  In  brief,  ZsGreen  retroviral 
particles  (Clontech)  were  isolated  from  PT-67  packaging  cells  (Clontech), 
filtered,  and  overlaid  onto  MCF-7  cells.  Cells  were  serially  transduced  twice 
for  24  h  each  round  and  then  subjected  to  G418  selection  followed  by 
aseptic  fluorescence-associated  cell  sorting  to  isolate  a  homogeneously 
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bright  green  subpopulation.  Generation  of  PR-B-expressing  T47D  human 
breast  cancer  cells  has  been  described  previously  (16). 

Xenograft  tumor  growth  and  metastases.  All  animal  procedures  were 
done  under  a  protocol  approved  by  the  University  of  Colorado  Institutional 
Animal  Care  and  Use  Committee.  Ovariectomized  female  athymic  nu/nu 
mice  were  obtained  from  Harlan  Sprague  Dawley  at  5  to  6  weeks  of  age. 
To  establish  tumors,  animals  were  anesthetized  with  Avertin  and  injected 
into  the  opening  of  the  lactiferous  duct  of  abdominal  mammary  glands  with 
1  million  ZsGreen-expressing  MCF-7  cells  in  100  pL  of  100%  Matrigel 
(Becton  Dickinson).  Mice  were  also  implanted  with  silastic  pellets  contain¬ 
ing  17(2)-estradiol  (2  mg  +  8  mg  cellulose)  as  previously  described  (17). 
Fluorescent  whole-body  imaging  (Illumatool  9900,  Lightools  Research)  was 
done  weekly  to  determine  when  axillary  LN  metastases  arose.  For  EWD 
mice,  once  axillary  LN  metastases  were  detected,  mice  were  anesthetized 
with  Avertin  and  estradiol-releasing  pellets  were  removed  for  1  week; 
some  were  removed  for  4  weeks  to  confirm  protein  expression  and  for 
proliferation  studies.  Control  mice  were  continued  on  estradiol.  Mice  were 
euthanized  by  C02  asphyxiation  and  fluorescent  intravital  optical  imaging 
was  done  by  coupling  the  Illumatool  with  an  Olympus  SZ-61  dissecting 
microscope  and  Olympus  C-5050  digital  camera.  Fluorescent  tumors  in 
mammary  glands  and  LNs  were  removed.  Lymphatic  vessels  (LV)  draining 
fluorescent  tumor  cells  from  mammary  glands  were  also  collected. 

Serial  passaging  of  a  ZsGreen+  T47D-PR-B  LN  metastasis  was  done  by 
isolating  the  LN,  pushing  the  tissue  through  multiple  sieves,  expanding  and 
G418  selecting  human  tumor  cells  in  culture,  and  reinjecting  the  cells  into 
the  mammary  glands  of  estradiol-supplemented  recipient  mice. 

Laser  capture  microdissection  and  expression  profiling.  Mammary 
gland  tumors  (MGT)  and  LN  metastases  from  estradiol-treated  or  1-week 
EWD  mice  were  harvested,  placed  in  a  Tissue-Tek  cryomold  (EMS),  covered 
with  Neg50  frozen  section  medium  (Richard- Allan  Scientific),  and  frozen 
in  isopentane  (Sigma)  cooled  by  liquid  nitrogen.  Sections  (8  pm)  were  cut 
through  the  entire  tissue  with  a  cryotome  set  at  —  20  °C.  Each  section  was 
collected  on  an  uncharged  slide  (Fisher),  and  every  10th  slide  was  stained 
for  H&E  or  processed  for  cytokeratin  18  (CK18)  immunofluorescence.  Slides 
containing  cancer  cells  were  processed  through  75%  ethanol,  water,  75% 
ethanol,  95%  ethanol,  and  100%  ethanol  for  30  s  each  and  finally  dehydrated 
with  xylene  for  1  min  before  being  immediately  subjected  to  laser  capture 
microdissection  using  an  Arcturus  Autopix.  Two  thousand  cells  for  each 
sample  were  extracted  within  30  min  of  fixation  and  frozen  at  —  70  °C  until 
further  processing.  Suggested  protocols  from  Arcturus  were  followed  for 
RNA  extraction  and  two  rounds  of  amplification  using  the  PicoPure  RNA 
isolation  kit  and  RiboAmp  HS  kit  (Arcturus),  respectively.  Samples  were 
in  vitro  transcribed,  biotinylated,  hybridized  to  Affymetrix  X3P  chips  that 
interrogate  47,000  transcripts  with  61,000  probe  sets,  and  scanned  on  an 
Agilent  Bioanalyzer  2100. 

Statistical  analyses.  Raw  expression  values  for  each  tumor  and  LN  were 
downloaded  into  GeneSpring  version  7.3  (Agilent)  for  exploratory  purposes. 
Four  matched  pairs  of  tumors  and  metastatic  LNs  (2  estradiol  and  2  EWD) 
as  well  as  2  unmatched  LNs  (1  estradiol  and  1  EWD)  for  a  total  of 
10  samples  were  analyzed.  Probe  sets  that  were  absent  across  all  10  samples 
were  filtered  from  further  statistical  analysis  (18).  Remaining  probe  sets 
were  robust  multiarray  analysis  (RMA)  normalized  in  Partek  Genomics 
Suite  6.2,  and  a  three-way  mixed  model  ANOVA  analysis  was  done  to  iden¬ 
tify  differentially  expressed  genes.  Four  groups  were  compared:  estradiol- 
treated  MGT,  estradiol-treated  LN,  EWD  MGT,  and  EWD  LN.  These  four 
groups  were  divided  based  on  the  two  levels  of  hormone  treatment 
(estradiol  or  EWD)  and  the  two  levels  of  tissue  type  (MGT  or  LN).  Hormone 
and  tissue  were  therefore  defined  as  main  factors  in  the  model.  Interac¬ 
tion  effects  among  hormone  and  tissue  were  also  included.  The  random 
factor  animal  (which  is  nested  in  hormone)  was  added  to  account  for  the 
fact  that  multiple  tissues  (MGT  and  LN)  were  taken  from  the  same  animal. 
In  total,  1,570  transcripts  were  identified  as  having  significant  differences 
among  the  four  group  means  (Pmodei  <  0.05)  and  P  <  0.05  for  main  or 
interaction  effects.  To  further  identify  differences  in  hormonal  regulation 
among  the  1,570  differentially  expressed  genes,  they  were  subjected  to  linear 
contrasts  to  make  two  specific  comparisons:  EWD  versus  estradiol-treated 
MGT  and  EWD  versus  estradiol-treated  LN.  Bonferroni’s  correction  for 


multiple  tests  was  then  applied  to  these  two  comparisons  and  the  level  of 
significance  was  changed  to  0.025. 

Gene  networking.  To  define  relationships  among  genes  differentially 
regulated  by  estradiol  in  MGTs  and  LN  metastases,  pathway  analyses  were 
done  using  Ingenuity  software.  Probe  set  identification  numbers  were  used 
for  gene  identification. 

Immunohistochemistry.  Organs  were  fixed  in  4%  paraformaldehyde 
overnight,  paraffin  embedded,  and  cut  into  4-pm  sections.  After  high- 
temperature  antigen  retrieval  in  citrate  buffer,  immunohistochemistries 
were  done  on  estradiol-treated  tissues  and  1-week  or  4-week  EWD  mice 
using  primary  antibodies  directed  against  ER  (1:100,  SP1;  Lab  Vision/ 
NeoMarkers),  progesterone  receptor  (PR;  1:500,  1294;  Dako),  CK18  (1:200; 
Calbiochem),  CD44  (1:200;  Ab4;  Lab  Vision/NeoMarkers),  caveolin  1  (CAV1; 
1:125;  Epitomics),  and  cathepsin  D  (CTSD;  1:125;  Epitomics)  for  1  h  at 
room  temperature.  Bound  primary  antibodies  were  detected  with 
horseradish  peroxidase  (HRP)-conjugated  goat  anti-mouse  and/or  goat 
anti-rabbit  secondary  antibodies  (Envision  HRP;  Dako)  reacted  with  3,3'- 
diaminobenzidine  (Dako).  Sections  were  counterstained  with  hematoxylin 
and  mounted  with  Permount  (Fisher).  For  immunofluorescence  microscopy, 
goat  anti-rabbit  Alexa  Fluor  555  (red)  secondary  antibodies  were  used 
(1:200,  Alexa  Fluors;  Invitrogen)  and  sections  were  counterstained  with 
4',6-diamidino-2-phenylindole  (DAPI)  in  methanol  and  mounted  with  Gel/ 
Mount  (Biomeda). 

Results 

EWD  model.  MCF-7+ZsGreen  tumors  require  estradiol  for 
growth  and  LN  metastasis  (15).  To  generate  LN  metastases, 
bilateral  tumors  were  grown  in  the  abdominal  mammary  glands 
of  ovariectomized  immunocompromised  mice  that  were  supple¬ 
mented  with  an  estradiol-releasing  pellet  (Fig.  L4).  Weekly  tumor 
measurements  were  recorded  with  a  digital  caliper.  Spread  of 
cancer  cells  to  axillary  LNs  was  determined  by  weekly  fluorescent 
whole-body  imaging  (Fig.  IB).  To  define  estradiol-regulated  genes 
(7),  once  LN  metastases  were  detected,  the  estradiol-releasing 
pellet  was  withdrawn  for  1  week  in  half  of  the  mice  (EWD)  or  left  in 
place  in  the  remaining  mice.  To  assess  proliferation  rate,  EWD  was 
extended  for  4  weeks  in  tumors  and  LN  metastases.  Bromodeox- 
yuridine  labeling  and  mitotic  indices  indicated  a  similar  decrease 
of  —50%  at  both  sites.  At  necropsy,  assessment  of  the  uterine 
mass  confirmed  the  presence  or  absence  of  systemic  estradiol,  with 
1-  and  4- week  EWD  uteri  exhibiting  significant  decreases  in  mass 
compared  with  estradiol-treated  uteri.  To  locate  cancer  cells  within 
LNs,  the  entire  LN  was  sectioned  and  every  10th  slide  was  pro¬ 
cessed  for  CK18  immunohistochemistry  (Fig.  1C).  The  CK18  map 
and  the  unique  morphology  of  MCF-7  cells  compared  with 
lymphocytes  (Fig.  ID)  allowed  for  identification  of  the  cancer  cells 
in  serial  sections.  We  then  marked  and  laser  captured  pure  MCF-7 
cell  populations  from  matched  tumors  and  LN  metastases. 

Estrogen  insensitivity  in  LN  metastases:  ER  and  PR 
expression.  Estradiol  treatment  leads  to  down-regulation  of  ER 
protein  levels  in  tumors  (19,  20).  To  characterize  these  effects  in 
estradiol-free,  estradiol-treated,  and  EWD  conditions,  the  MCF-7 
tumors,  LV  emboli,  and  LN  metastases  were  probed  for  ER  and 
PR  by  immunohistochemistry.  In  Fig.  2,  the  tissues  from  each 
treatment  group  were  from  the  same  mouse.  Highest  ER  levels 
were  present  in  tumors  that  did  not  receive  supplemental  estradiol 
(— E).  Seventy- three  percent  of  MCF-7  cells  expressed  high  levels 
of  ER  and  no  cells  contained  PR  (Fig.  2,  top).  In  the  presence 
of  estradiol  (+E),  ERs  were  down-regulated  by  >40%  in  the  tumor 
(Fig.  2,  middle)  but  decreased  by  <20%  in  LV  tumor  emboli  and  LN 
metastasis  of  the  same  mouse.  Ligand-dependent  down-regulation 
of  ER  is  required  for  its  transcriptional  activity  (20),  including  PR 
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Figure  1.  MCF-7+ZsGreen  xenograft  tumors  and  LN  metastases. 
MCF-7+ZsGreen  tumors  were  established  in  mammary  glands  of  ovariectomized 
nu/nu  mice  supplemented  with  an  estradiol  pellet,  and  LN  metastases 
were  allowed  to  develop.  A,  fluorescent  whole-body  image  of  bilateral 
MCF-7+ZsGreen  tumors.  B,  fluorescent  whole-body  image  of  MCF-7+ZsGreen 
metastasis  to  the  left  axillary  LN.  C,  fluorescence  immunohistochemistry  for 
CK18  {red)  and  DAPI  {blue)  of  an  axillary  LN  metastasis.  Magnification,  x4. 

D,  white  light  laser  capture  microscope  image  of  an  axillary  LN  metastasis 
showing  MCF-7  cells  or  lymphocytes  (LN).  Magnification,  x40. 


induction.  Indeed,  PRs  were  expressed  in  +E  tumors  but  were  only 
rarely  found  in  +E  LV  emboli  (Fig.  2,  arrows )  or  in  LN  metastases. 
On  EWE),  ERs  were  restored  to  control  (— E)  levels  in  all  tissues,  and 
PRs  were  absent  (Fig.  2,  bottom ). 


Microenvironment  versus  clonal  selection.  To  determine 
whether  the  increased  ER  expression  observed  in  LN  metastases 
was  due  to  clonal  selection  of  ER+  cells  from  the  primary  tumor,  or 
to  an  effect  on  ER  at  the  metastatic  microenvironment,  tumor  cells 
were  isolated  from  the  axillary  LN  metastasis  of  a  donor  mouse, 
expanded  in  culture,  and  reimplanted  in  the  mammary  gland  of  an 


Figure  2.  Estradiol  regulation  of  ER  and  PR  levels.  ER  and  PR  expression 
in  MCF-7  cells  from  matched  tumors,  LVs,  and  LNs.  MCF-7+ZsGreen  tumors 
were  established  in  ovariectomized  mice  supplemented  with  a  placebo  pellet 
(-E),  an  estradiol-releasing  pellet  (+E),  or  an  estradiol-releasing  pellet  until  LN 
metastases  arose,  after  which  the  estradiol  pellet  was  withdrawn  for  1  wk 
(El/I/D).  The  primary  tumors,  efferent  LVs,  and  draining  LNs  from  the  same 
mouse  were  dissected  and  prepared  for  histology.  Paraffin-embedded  sections 
were  probed  with  rabbit  anti-ER  or  mouse  anti-PR  primary  antibodies  coupled 
with  goat  anti-rabbit  (ER)  or  goat  anti-mouse  (PR)  HRP-conjugated  secondary 
antibodies.  Background  staining  of  erythrocytes  (e).  Arrows  point  to  the  LV. 
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Figure  3.  Aberrant  ER  regulation  in  tumor 
cells  located  in  the  LN  microenvironment 
and  global  gene  expression  profiling 
of  tumors  and  matched  LN  metastases. 

A,  ER  immunohistochemistry  of  estradiol- 
treated  T47D  tumors  and  their  matched  LN 
metastases.  An  estradiol-treated  MGT  in 
mouse  1  developed  LN  metastases. 

Cells  from  the  LN  metastasis  of  mouse  1 
were  expanded  in  culture  and  reinjected 
into  the  mammary  gland  of  a  second 
estradiol-treated  recipient  mouse  2, 
which  also  yielded  LN  metastases.  MGTs 
and  LN  metastases  from  both  mice 
were  paraffin  embedded,  sectioned,  and 
stained  for  ER  by  immunohistochemistry. 

B,  unsupervised  hierarchical  clustering 
from  continuously  estradiol-treated  (+E) 
or  1-wk  EWD  tumors  and  matched 

LN  metastases.  Red,  up-regulated 
transcripts;  blue,  down-regulated 
transcripts;  gray,  relatively  unchanged 
transcripts.  C,  ANOVA  analysis  identified 
1 ,570  genes  with  tissue  effects,  hormone 
effects,  and/or  interaction  effects. 
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estradiol-treated  recipient  mouse.  This  tumor  also  yielded  a  LN 
metastasis.  Figure  3A  {top)  shows  ER  immunohistochemistry  of 
the  estradiol-treated  primary  MGT  and  its  matched  LN  metastasis 
from  the  donor  mouse  (mouse  1).  ERs  in  the  tumor  grown  in  the 
mammary  gland  and  LN  of  the  recipient  mouse  (mouse  2)  are 
shown  in  Fig.  3A  {bottom).  There  was  a  2-fold  increase  in  ER 
content  in  the  cancer  cells  populating  the  LN  of  mouse  1  and  a 
similar  2-fold  increase  in  the  LN  ER  of  mouse  2.  These  data  suggest 
that  the  LN  microenvironment  rather  than  clonal  selection 
modifies  ER  expression  levels  in  response  to  estradiol. 

Global  gene  expression  profiling.  To  determine  if  the  apparent 
estradiol  insensitivity  of  LN  metastases  as  assessed  by  ER  down- 
regulation  and  failure  of  PR  induction  extends  to  other  estradiol- 
dependent  genes,  global  gene  expression  profiling  was  done.  To 
ensure  that  only  human  tumor  cells  were  evaluated,  pure  cell  popu¬ 
lations  were  isolated  by  laser  capture  microdissection.  Estradiol- 
treated  and  1-week  EWD  tumors  and  their  matched  LN  metastases 
were  analyzed  to  define  the  subset  of  estradiol-regulated  genes  (7). 
Unsupervised  hierarchical  clustering  showed  that  EWD  had  a  lesser 
effect  on  genes  in  LN  metastases  than  on  genes  in  MGTs,  with  the 
estradiol-treated  and  EWD  LN  metastases  clustering  closest 


together  (Fig.  3 B).  Interestingly,  the  EWD  tumors  clustered  closer 
to  the  LNs  than  to  the  estradiol-treated  tumors. 

To  identify  genes  with  differential  expression  patterns,  a  three- 
way  mixed  model  ANOVA  was  done  for  genes  that  were  present 
in  at  least  one  condition.  Genes  with  the  greatest  differences  in 
expression  due  to  the  hormone  treatment,  tissue  type,  interaction 
between  hormone  and  tissue,  or  within  and  between  animals 
(Pmodei  <  0*05)  and  with  significant  (P  <  0.05)  main  (tissue  and/ 
or  hormone  effects)  or  interaction  effects  were  sorted  to  yield  a 
total  of  1,570  differentially  expressed  transcripts.  Significant  main 
effects  indicate  simple  differences  between  the  two  hormone  and/ 
or  two  tissue  groups,  an  example  being  that  all  EWD  tissues  have 
higher  gene  expression  levels  than  all  estradiol-treated  tissues 
(regardless  of  tissue  type)  or  that  all  LN  tissues  have  higher  gene 
expression  levels  than  all  MGT  (regardless  of  hormone  treatment). 
Significant  interaction  effects,  on  the  other  hand,  indicate  com¬ 
plex  differences  in  which  the  results  of  one  factor  differ  depending 
on  the  levels  of  the  other  factor.  For  example,  a  gene  that  was 
expressed  at  higher  levels  in  estradiol-treated  LNs  than  MGTs  but 
in  EWD  tissues  had  the  opposite  pattern  with  higher  levels  in 
MGTs  than  LNs. 
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Figure  4.  Genes  that  are  tissue  and 
hormone  regulated  with  or  without 
interaction  effects.  RMA  normalized  mean 
RNA  expression  values  with  SE.  Solid 
lines,  data  from  estradiol-treated  tissues; 
dashed  lines,  data  from  EWD  tissues. 

A,  tissue-regulated  genes.  B,  hormone- 
regulated  genes.  MT1F,  metallothionein 
1 F.  C,  genes  that  are  tissue  and  hormone 
regulated.  D,  genes  with  interaction  effects 
that  are  differentially  regulated  by  estradiol 
in  the  two  tissues  include  CAV1,  MUC5AC, 
TP53,  BCL2,  CDKN1A,  and  CD44. 


Of  the  1,570  differentially  expressed  transcripts,  more  had 
increased  average  expression  in  LN  metastases  (70%)  than  in 
primary  tumors  (30%).  In  addition,  on  average,  genes  were  more 
highly  expressed  in  EWD  tissues  (64%)  than  in  estradiol-treated 
ones  (36%),  suggesting  that  estradiol  affects  approximately  two 
thirds  of  genes  through  suppression  of  RNA  levels.  Among  the 
1,570  transcripts,  386  had  only  significant  (P  <  0.05;  Fig.  3C)  tissue 
effects,  44  had  only  significant  hormone  effects,  and  91  were 
affected  by  both  the  tissue  microenvironment  and  the  presence  of 
hormones  (Fig.  3 C).  However,  the  majority  of  genes,  1,049  of  1,570, 
exhibited  a  significant  interaction  effect.  The  expression  values 
of  these  genes  could  not  be  simply  described  as  due  to  either  of 
the  main  effects  but  rather  due  to  an  interrelationship  between 
them.  Most  of  these  genes  may  also  have  significant  main  effects  as 
shown  by  intersection  in  the  Venn  diagram,  but  interaction  effects 
take  precedence  over  these  main  effects,  and  the  main  effects 
cannot  be  correctly  interpreted  without  acknowledging  the  inter¬ 
action  effects  (21).  For  these  genes,  expression  levels  in  the  two 
tissues  (MGT  and  LN)  were  influenced  differently  by  the  type  of 


hormone  treatment  (estradiol  and  EWD).  A  subset  (145  of  1,049)  of 
these  genes  exhibited  only  interaction  effects. 

Main  effects  and  interaction  effects  can  be  visualized  by 
graphical  representation  with  levels  of  one  factor  on  the  X  axis 
and  separate  lines  representing  levels  of  the  other  factor,  as  in 
Fig.  4.  Lines  that  are  parallel  to  one  another  with  steep  slopes  or 
separation  indicate  main  effects,  whereas  nonparallel  lines  indicate 
interaction  effects.  Figure  4 A  shows  two  examples  of  genes  with 
significant  tissue  effects.  Cortactin  transcripts  had  significantly 
higher  expression  levels  in  LN  metastases,  regardless  of  hormone, 
as  did  the  majority  (290  of  386)  of  the  tissue-regulated  genes  (other 
examples  include  MAPK  and  cyclin  F).  CXCL14  is  one  of  the  96  of 
386  tissue-regulated  genes  whose  levels  were  higher  in  the  primary 
tumors  compared  with  LN  metastases. 

Figure  4 B  shows  two  examples  of  genes  with  a  significant  hor¬ 
mone  effect.  The  expression  of  cathepsin-D  (CTSD)  transcripts  was 
dependent  on  estradiol  regardless  of  the  microenvironment.  Other 
genes  with  similar  expression  patterns  were  insulin-like  growth 
factor-I  receptor  ( IGF-IR ),  cadherin-15,  and  defensin  [Fl.  Transcript 
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levels  of  metallothionein  IF  were  suppressed  by  estradiol  regardless 
of  the  microenvironment.  Other  examples  of  estradiol-suppressed 
genes  are  cathepsin  F  and  BLNK. 

A  subset  of  differentially  expressed  genes  (91  of  1,570)  had  RNA 
expression  levels  that  were  both  hormone  and  tissue  regulated  as 
shown  in  Fig.  4C.  ERBB3  was  significantly  suppressed  by  estradiol 
regardless  of  the  microenvironment,  but  its  expression  was  also 
increased  in  LN  metastases  compared  with  MGTs  of  both  hormone 
treatments.  SEBP1  was  also  significantly  increased  in  LN 
metastases  compared  with  MGTs  regardless  of  the  presence  of 
estradiol  and  was  also  significantly  up-regulated  by  estradiol  in 
both  tissues. 

The  majority  of  differentially  expressed  genes  (1,049  of  1,570)  had 
a  significant  interaction  effect,  exhibiting  an  interrelationship 
between  tissue  and  hormonal  factors.  Figure  4 D  shows  examples  of 
a  few  such  genes.  Transcripts  for  BCL2,  CAV1,  CD44,  CDKN1A, 
MUC5AC,  and  tumor  promoter  53  ( TP53 )  all  exhibit  interaction 
effects  and  are  affected  differently  by  estradiol  in  tumors  and  LN 
metastases.  MUC5AC  and  TP53  are  examples  of  genes  exhibiting 
main  effects  (hormone  and  tissue)  but  with  the  interaction  term 
predominating.  Main  effects  of  CDKN1A  (p21 )  cancel  each  other 
out,  leaving  the  dominant  role  of  interaction  effects  evident. 

Immunohistochemistry.  Several  of  the  genes  that  were 
identified  by  the  expression  profiling  analysis  were  investigated 
further  at  the  protein  level  by  immunohistochemistry  (Fig.  5).  The 
estradiol-regulated  gene  CTSD  was  expressed  in  80%  to  90%  of 
estradiol- treated  MCF-7  primary  tumor  cells,  and  EWD  reduced 
this  number  to  50%  to  60%.  In  the  LN  metastases  in  the  presence  of 
estradiol,  50%  to  60%  of  cancer  cells  expressed  CTSD;  this  fell  to 
10%  to  15%  on  EWD  (Fig.  5,  top).  CAV1  protein  was  unaffected  by 
EWD  but  was  expressed  only  in  MGTs  and  not  in  LN  metastases 
(Fig.  5,  middle ).  In  estradiol-treated  tumors,  CD44  was  heteroge¬ 
neously  expressed;  on  EWD,  nonsignificant  increases  in  the 
percentage  of  CD44+  cells  were  recorded.  Interestingly,  the  opposite 
pattern  was  observed  in  LN  metastases.  The  extensive  over¬ 
expression  of  CD44  in  estradiol- treated  cancer  cells  metastatic  to 
LNs  (15)  is  dramatically  reduced  by  EWD  (Fig.  5,  bottom ). 

Estrogen-regulated  genes.  Two  linear  contrasts  were  done 
among  the  1,570  differentially  expressed  transcripts  comparing 
estradiol-treated  with  EWD  tumors  and  estradiol-treated  with 
EWD  LN  metastases  (Fig.  6 A  and  B).  This  identified  transcripts  that 
were  estradiol  regulated  in  MGTs  compared  with  LN  metastases. 
There  were  273  significantly  hormonally  up-regulated  or  down- 
regulated  genes  (fold  change  >1.3;  P  <  0.025;  Fig.  6 A).  Of  these,  256 
were  estradiol  dependent  in  tumors,  41  were  estradiol  dependent  in 
LN  metastases,  and  24  maintained  estradiol  sensitivity  in  both 
MGTs  and  LN  metastases  (Fig.  6 A).  These  data  confirm  the  con¬ 
clusions  drawn  from  Figs.  2  and  3 A  that  tumor  cells  in  the  LN 
microenvironment  are  relatively  estradiol  resistant. 

To  identify  genes  with  any  absolute  difference  in  expression  due 
to  estradiol  in  MGTs  and  LN  metastases,  the  fold  change  criterion 
was  removed  from  the  two  linear  contrasts.  This  increased  the 
total  list  of  significantly  (P  <  0.025)  hormone-regulated  genes  from 
273  to  1,176  (Fig.  6 B).  These  additional  903  genes  were  up-regulated 
or  down-regulated  with  fold  changes  >1.0  and  <1.3,  and  in  total, 
there  were  707  estradiol-regulated  genes  in  MGTs,  177  estradiol- 
regulated  in  LN  metastases,  and  292  estradiol-regulated  in  both 
sites.  Seventy-five  of  these  transcripts  were  regulated  by  estradiol 
in  opposite  directions  in  MGTs  and  LN  metastases  (Fig.  6C). 
BCL2,  CD44,  and  CDKN1A  are  examples  of  such  genes  (Fig.  4 D). 
CDKN1A,  also  known  as  p21/WAF,  encodes  a  potent  cyclin- 
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Figure  5.  Immunohistochemical  analyses  of  proteins  that  are  hormone  and 
tissue  regulated.  Estradiol-treated  and  EWD  tumors  and  matched  LN 
metastases  from  the  same  mouse  were  paraffin  embedded,  sectioned,  and 
probed  with  CTSD,  CAV1,  or  CD44  primary  antibodies  coupled  with  goat 
anti-rabbit  or  goat  anti-mouse  HRP-conjugated  secondary  antibodies.  Top, 
estradiol-regulated  CTSD;  middle,  CAV1  protein  regulation  by  the  tissue 
microenvironment;  bottom,  differential  estradiol  regulation  of  CD44  expression 
in  tumors  and  LN  metastases. 

dependent  kinase  inhibitor  and  regulator  of  cell  cycle  progression 
at  Gx.  Its  expression  is  tightly  controlled  by  TP53,  one  of  our  set 
of  microenvironment-regulated  genes  (Fig.  4 D).  CD44  (Fig.  5)  is  a 
putative  breast  cancer  stem  cell  marker  (22)  and  may  target  tumor 
cells  to  LNs  (15). 

Network  analysis.  Ingenuity  data  analysis  was  conducted  to 
determine  if  the  genes  that  had  interaction  effects  and  were 
differentially  regulated  by  estradiol  in  MGTs  and  LN  metastases  are 
also  involved  in  important  functional  networks.  Twenty-seven 
unique  gene  networks  that  contained  at  least  one  differentially 
regulated  gene  were  identified.  Supplementary  Fig.  SI  shows  a 
network  of  33  genes,  100%  of  which  are  differentially  estradiol 
regulated  in  MGTs  and  LN  metastases.  This  network  has  TP53  at 
its  center  and  functions  in  cellular  responses  to  “therapeutics, 
cell  death,  and  connective  tissue  disorders.”  Most  of  these  genes 
(25  of  33)  were  estradiol  regulated  in  MGTs  but  lost  estradiol  sen¬ 
sitivity  in  LN  metastases  (orange  outline).  Two  genes  were  estradiol 
regulated  in  LN  metastases  but  not  in  tumors  (red  outline),  and  the 
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remaining  seven  genes  were  regulated  in  the  opposite  direction  by 
estradiol  in  MGTs  and  LN  metastases  (blue  outline). 

Discussion 

Gene  expression  changes  in  MGTs  and  LN  metastases.  Two 

major  theories  have  been  advanced  to  explain  successful  cancer 
metastasis,  (a)  The  “rare  tumor  cell  hypothesis”  states  that  only 
certain  cells  within  tumors  are  capable  of  migrating  and/or 
growing  at  distant  sites.  These  rare  cells  may  be  stem  cells,  such  as 
those  that  express  CD44  (22).  On  the  other  hand,  they  may 
represent  rare  cells  within  a  population  of  genetically  diverse  and 
unstable  cells  that  originated  with  or  acquired  the  appropriate 
genetic  information  needed  to  seed  at  specific  metastatic  sites  (23). 
(b)  An  alternative  hypothesis  posits  that  metastatic  capacity  is  an 
inherent  rather  than  an  acquired  feature  of  breast  cancers  (24,  25) 
and  is  based  on  data  suggesting  that  metastatic  capacity  is 
encoded  early  in  the  majority  of  cells  within  a  tumor.  In  support 
of  this  idea  are  studies  showing  that  tumors  and  their  LN  meta¬ 
stases  retain  similar  genetic  signatures,  including  aggressiveness 
and  proliferation  markers  (24,  26,  27).  One  study,  for  example, 
found  only  27  changed  genes  between  lung  tumors  and  their  LN 
metastases  (28).  Weigelt  et  al.  (27)  identified  a  70-gene  cluster  in 
primary  breast  tumors  that  predicts  likelihood  of  subsequent 
metastases  because  the  unique  expression  profiles  of  these  genes 
were  maintained  in  LNs  and  more  distant  metastases.  Another 
study  found  no  “shared  genes”  between  primary  tumors  that  were 
predictive  of  LN  metastases  but,  anywhere  between  3  and  149 
genes,  were  “antiexpressed”  or  expressed  in  opposite  directions  in 
tumors  versus  LN  metastases  (29). 

Our  data  provide  evidence  in  support  of  both  theories.  In  these 
studies,  we  analyze  the  changes  in  gene  expression  induced  by 
estradiol  as  revealed  by  estradiol  removal.  Fewer  than  5%  of  all 
genes  analyzed  were  differentially  expressed  in  cancer  cells  within 
the  mammary  gland  compared  with  the  LN.  Although  this  might  be 
considered  a  subtle  difference,  it  nevertheless  represents  differen¬ 
tial  regulation  of  >1,000  genes  in  these  two  microenvironments. 
These  genetic  differences  could  be  sufficient  to  define  a  “metastatic 
stem  cell”  (30),  expressing  the  adaptive  proteins  required  for 
successful  colonization  and  growth  of  a  subpopulation  in  a  meta¬ 
static  microenvironment,  as  envisioned  by  the  rare  tumor  cell 
hypothesis.  Elevated  expression  levels  of  the  hyaluronan  receptor 
CD44,  a  putative  cancer  stem  cell  marker,  in  estradiol-treated  LN 
metastases  suggest  such  selectivity  in  our  models.  On  the  other 
hand,  our  observations  that  cancer  cell  emboli  move  in  bulk 
through  LVs  to  LNs  (15),  coupled  with  our  findings  here  and  those 
of  others  (24,  27,  29)  of  modest  genetic  changes  between  primary 
tumors  and  LN  metastases,  could  be  interpreted  as  supporting 
the  view  that  the  majority  of  primary  tumor  cells  arise  either  with 
or  without  the  predisposition  for  LN  spread.  If  true,  the  gene 
expression  changes  we  observe  could  be  attributed  to  the  plasticity 
that  all  tumor  cells  exhibit  when  growing  in  two  distinct  micro¬ 
environments.  Evidence  supporting  these  microenvironmental 
effects  on  ER  expression  and  function  was  found  when  a  down- 
regulation  “resistant”  ER+  LN  metastasis  regained  ER  down- 
regulation  capacity  when  grown  in  the  mammary  gland  of  a 
recipient  mouse  (Fig.  3 A). 

ER  function  in  metastases.  Beatson  (31)  reported  >100  years 
ago  that  an  essential  feature  of  breast  cancer  is  continuous  and 
excessive  growth  of  the  epithelium,  which  invades  the  surround¬ 
ing  tissues,  spreads  along  the  LVs,  passes  from  one  set  of  LNs  to 


another,  and  eventually  forms  deposits  in  distant  organs,  steps  that 
occur  more  rapidly  and  prove  to  be  more  quickly  fatal  in  younger 
patients.  Knowing  that  the  ovaries  produced  factors  required  for 
lactation,  and  believing  “that  all  pathological  changes  are  merely 
modified  physiological  ones,”  Beatson  was  the  first  to  find  that 
removal  of  the  ovaries  and  the  estradiol  they  elaborate  reduced 
breast  tumor  size  and  extended  the  survival  of  young  women  with 
advanced  breast  cancer. 

Seventy  percent  to  80%  of  primary  breast  cancers  are  ER+,  and 
—  80%  of  LN  metastases  retain  the  receptors  (11-13,  32).  If  such 
tumors  are  treated  with  antiestrogens  or  EWD  therapies,  20%  to 
50%  acquire  resistance  (33),  often  recurring  as  LN  metastases. 
Premenopausal  women  with  LN  metastases  are  less  responsive  to 
tamoxifen  treatment  than  LN-negative  patients  (34),  and  nodal 
involvement  is  a  significant  predictor  of  local  relapse  (35)  with  a 
5.6-fold  higher  relative  risk  of  local  recurrence  than  women  with 
clear  LNs.  Preoperative  chemotherapy  causes  significantly  less 
damage  to  LN  metastases  in  comparison  with  primary  tumors  (36). 
Based  on  these  data,  we  hypothesized  that  ER  effects  on  gene 
expression  would  be  altered  by  the  LN  microenvironment.  To  the 
best  of  our  knowledge,  no  previous  studies  have  defined  estradiol- 
regulated  genes  in  tumors  compared  with  their  matched  LN  meta¬ 
stases.  Identifying  these  genes  is  an  important  step  toward 
elucidating  mechanisms  that  contribute  to  breast  cancer  progres¬ 
sion  and  estradiol  resistance  in  advanced  disease. 

In  the  face  of  continuous  estradiol  treatment  in  two  separate  ER+ 
breast  cancer  cell  lines,  higher  levels  of  ER  were  observed  in  LN 
metastases  than  in  matched  MGTs  (Figs.  2  and  3 A).  The  partial 
failure  of  ligand-dependent  ER  down-regulation  in  LN  metastases 
suggests  that  the  receptors  are  relatively  estradiol  insensitive  in  the 
LN  microenvironment.  Reduced  PR  expression  in  estradiol-treated 
LN  metastases  supports  this  conclusion,  as  does  global  gene 
expression  profiling,  which  showed  that  estradiol  had  less  influence 
on  gene  expression  in  LN  metastases  than  in  matched  primary 
MGTs  (Figs.  3  and  6).  We  have  now  identified  numerous  genes  that 
are  estradiol  regulated  in  primary  tumors  but  not  in  LN  metastases. 
Examples  include  MUC5A  and  TP53,  which  were  estradiol  sensitive 
in  MGTs  but  lost  this  regulation  in  the  LNs.  Additionally,  75  genes 
were  regulated  in  the  opposite  direction  by  estradiol  in  MGTs 
compared  with  their  LN  metastases  (Fig.  6C),  21  of  which, 
including  BCL2  and  CDKN1A  (Fig.  4),  were  previously  reported  to 
be  estradiol  regulated  (7,  37,  38).  CD44  (Figs.  4  and  5)  was  also 
regulated  in  opposite  directions  at  the  two  sites.  It  is  one  of  two  cell 
surface  markers  that  differentiate  tumorigenic  from  nontumori- 
genic  breast  cancer  cells  (22).  This  putative  stem  cell  and  aggres¬ 
siveness  marker  (39,  40)  is  an  example  of  how  EWD  differentially 
alters  the  expression  of  a  potential  therapeutic  target  in  primary 
tumors  compared  with  LN  metastases. 

Mechanisms  of  estradiol  resistance  in  LNs.  Possible  explan¬ 
ations  for  the  differential  estradiol  sensitivity  of  MGTs  versus  LN 
metastases  include  increased  availability  of  growth  factors  in  LNs 
that  directly  interfere  with  ER  transcription  (41),  diminished 
ubiquitination  and  loss  of  proteasome-mediated  protein  down- 
regulation  in  LNs  (42,  43),  loss  of  CAV1  and  subsequent  membrane 
signaling  (44),  and  reduced  estradiol  availability  to  LNs.  With 
regard  to  growth  factors,  IGFs  and  epidermal  growth  factors  (EGF) 
activate  the  phosphatidylinositol  3-kinase-AKT-mammalian  target 
of  rapamycin  pathways,  which  alter  Jun  and  Fos  activities  at  the 
PR  promoter  and  decrease  its  ER-mediated  transcription  (45). 
LN  stromal  cells  produce  both  IGF  and  EGF  (41),  which  may 
explain  the  reduced  PR  induction  as  well  as  the  reduced  estradiol 
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1 
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0.003 
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PIGH 
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D 

1 
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C12orf22 
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D 

1 
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POLD4 

0.011 

1 

D 

AU 148054 

C14orf82 

0.002 

1 

D 

NM_007198 

PROSC 

0.014 

1 

D 

NM_0 18840 

C20orf24 

0.024 

D 

1 

NM_002802 

PSMC1 

0.012 

D 

1 

NM_006429 

CCT7 

0.015 

D 

1 

NM_1 52442 

RAD9B 
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D 

1 

BC004372 

CD44 

0.002 

D 

1 

NM_0 15523 

REX02 

0.008 

D 

1 

U07969 

CDH17 

0.004 

D 

1 
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RGS3 
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1 

D 

NM_000389 

CDKN1A 

0.005 

D 

1 

AI884906 

RNF182 

0.005 

D 

1 

BE887449 

CROP 

0.016 

1 

D 

NM_000990 

RPL27A 

0.012 

D 

1 

NM  001896 

CSNK2A2 

0.012 

D 

1 

AL1 18510 

RPS10 
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1 
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DDOST 
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D 

1 

BC004886 

RPS17 

0.008 

D 

1 

NM_004090 

DUSP3 

0.011 

1 

D 
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RPS24 

0.006 

D 

1 

BC002719 

EIF3S1 
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D 

1 
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1 

D 

BG026723 
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D 

1 

AL1 36551 
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1 
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FLJ33996 
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1 
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D 

1 
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1 
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0.011 
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1 
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1 
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1 
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1 
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1 
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0.007 

1 

D 
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TBPIP 
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D 

1 

AU 1441 36 

KIFAP3 

0.003 

1 

D 

NM_004608 

TBX6 

0.007 

1 

D 

L38951 

KPNB1 

0.010 

D 

1 

AL1 33033 

THRAP2 

0.013 

D 

1 

NM  000421 

KRT10 

0.005 

D 

1 

AC007954 

TMED8 

0.001 

1 

D 

NM_002289 

LALBA 

0.010 

1 

D 

AA470128 

TYSND1 

0.021 

1 

D 

AU 145563 

LOC255326 

0.005 

1 

D 

NM_024102 

WDR77 

0.012 

D 

1 

BC043252 

LOC283874 

0.009 

1 

D 

NM_0 16384 

ZMYM5 

0.018 

D 

1 

BF683512 

LOC284184 

0.001 

D 

1 

AW014719 

ZNF343 

0.003 

1 

D 

BE551864 

LOC440149 

0.005 

1 

D 

Figure  6.  Pairwise  comparisons  to  identify  estradiol-regulated  transcripts  in  tumors  (T),  LN  metastases,  or  both.  A,  transcripts  regulated  in  a  significant  manner 
(P  <  0.025)  by  estradiol  >1. 3-fold.  B,  transcripts  regulated  in  a  significant  manner  (P  <  0.025)  by  estradiol  regardless  of  fold  levels.  C,  list  of  75  transcripts  that  are 
regulated  in  a  significant  manner  in  opposite  directions  by  estradiol  in  MGTs  compared  with  LN  metastases.  /,  increases;  D,  decreases  with  estrogen. 


regulation  of  other  target  genes  at  this  site.  Increased  cortactin 
levels  (Fig.  4)  inhibit  ubiquitin-mediated  degradation  of  EGF 
receptors  following  their  internalization,  resulting  in  sustained 
ligand-induced  EGF  signaling  (42,  43).  It  is  possible  that  similar 


ubiquitin-dependent  mechanisms  explain  the  ER  down-regulation 
defects  we  observe  in  LNs  (Fig.  2;  ref.  15).  Supporting  this  hypo¬ 
thesis,  nine  transcripts  known  to  be  involved  in  the  ubiquitination 
process  were  differentially  expressed  in  MGTs  and  LN  metastases. 
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This  includes  two  known  enzymes  that  function  in  the  ubiquitina- 
tion  process,  UBE2D3  and  USP7,  which  were  increased  in  LN 
metastases  and  also  lost  estradiol  sensitivity  in  LN  metastases 
(Supplementary  Fig.  SI).  CAV1,  a  scaffold  protein  that  transduces 
ER  signaling  at  the  plasma  membrane  (44),  was  decreased  in  LNs, 
possibly  explaining  loss  of  ER  responsiveness.  Finally,  although  it  is 
possible  that  the  diminished  estradiol  sensitivity  of  LN  meta¬ 
stases  is  due  to  an  overall  decrease  of  estradiol  levels  within  LNs, 
this  is  highly  unlikely  given  the  number  of  estradiol-regulated  genes 
whose  levels  are  maintained  or  increased  in  LN  metastases  (Fig.  4). 

Networks.  We  speculated  that  genes  differentially  regulated  in 
MGTs  versus  LN  metastases  may  contribute  to  tumor  aggressive¬ 
ness,  recurrence,  and/or  growth  in  LNs.  We  therefore  analyzed  the 
functional  networks  characterizing  these  genes.  Of  interest  was  a 
33-gene  network  composed  entirely  of  transcripts  regulated 
differently  by  estradiol  in  MGTs  and  LNs  (Supplementary  Fig.  SI), 
which  included  CDNK1A,  TP53,  and  BCL2-linked  signaling  path¬ 
ways.  In  this  network,  the  majority  of  the  genes  were  significantly 
regulated  by  estradiol  in  tumors  while  losing  estradiol  sensitivity  in 
LN  metastases  or  suppressed  by  estradiol  in  tumors  but  increased 
by  estradiol  in  LN  metastases.  These  findings  are  interesting  given 
the  fact  that  70%  of  the  1,570  global  genes  differentially  expressed 
in  the  two  sites  were  increased  in  LN  metastases  compared  with 
MGTs,  whereas  64%  of  the  differentially  estradiol-regulated  genes 
were  suppressed  by  estradiol.  This  corroborates  previous  data  with 
MCF-7  cells  showing  that  70%  of  genes  are  down-regulated  by 
estradiol  treatment  (38)  and  is  consistent  with  an  overall  lifting  of 
the  suppressive  effects  of  estradiol  in  MCF-7  tumor  cells  growing  in 
LNs.  It  is  possible  that  these  results  do  not  completely  mimic 
effects  in  a  patient  with  an  intact  immune  system.  The  LNs  of  nude 
mice  lack  or  have  reduced  numbers  of  T  cells  (46).  However,  they 
have  B  cells,  natural  killer  cells,  reticular  cells,  and  macrophages  as 
well  as  a  nodal  stromal  network  and  architecture.  Because  different 


substrata  are  known  to  alter  growth  and  function  of  ER+  human 
breast  cancer  cells  in  vitro  (6),  it  is  likely  that  the  LN  micro¬ 
environment,  including  its  unique  architecture,  cytokines,  and 
growth  factors,  modifies  ER  function.  Given  the  experimental 
nature  of  our  models,  these  findings  stress  the  need  for  further 
studies  to  show  whether  estradiol  insensitivity  occurs  in  clinical 
ER+  LN  metastases  and  whether  other  metastatic  microenviron¬ 
ments  differentially  modulate  estradiol  sensitivity  and  to  define 
estradiol-regulated  genes  uniquely  influenced  by  the  microenvi¬ 
ronment  of  tumor  cells.  This  approach  may  identify  key  genes  that 
can  be  targeted  therapeutically  at  specific  metastatic  sites. 

In  summary,  we  describe  a  xenograft  model  of  EWD  after  LN 
metastasis  and  show  that  estradiol  regulation  of  genes  is  different 
in  ER+  primary  MGTs  compared  with  their  matched  ER+  LN 
metastases.  Through  global  gene  expression  profiling  and  immu- 
nohistochemical  analyses,  we  identify  genes  that  maintain,  lose,  or 
are  regulated  in  the  opposite  direction  by  estradiol  in  tumors  and 
LN  metastases.  These  results  suggest  that  current  ER-targeted 
therapies  may  differentially  affect  tumor  cells  depending  on  the 
microenvironment  in  which  the  cells  reside.  If  so,  our  studies 
provide  possible  new  therapeutic  targets,  which,  in  combination 
with  standard  hormonal  therapies,  might  improve  responses  of 
advanced  breast  cancers. 
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Abstract  Genome-wide  expression  profiling  has  expe¬ 
dited  our  molecular  understanding  of  the  different  subtypes 
of  breast  cancers,  as  well  as  defined  the  differences  among 
genes  expressed  in  primary  tumors  and  their  metastases. 
Laser-capture  microdissection  (LCM)  coupled  to  gene 
expression  analysis  allows  us  to  understand  how  specific 
cell  types  contribute  to  the  total  cancer  gene  expression 
signature.  Expression  profiling  was  used  to  define  genes  that 
contribute  to  breast  cancer  spread  into  and/or  growth  within 
draining  lymph  nodes  (LN).  Whole  tumor  xenografts  and 
their  matched  whole  LN  metastases  were  compared  to  LCM 
captured  cancer  cells  from  the  same  tumors  and  matched 
LN  metastases.  One-thousand  nine-hundred  thirty  genes 
were  identified  by  the  whole  organ  method  alone,  and  1,281 
genes  by  the  LCM  method  alone.  However,  less  than  1% 
(30  genes)  of  genes  that  changed  between  tumors  and  LN 
metastases  were  common  to  both  methods.  Several  of  these 
genes  have  previously  been  implicated  in  cancer  aggres¬ 
siveness.  Our  data  show  that  whole-organ  and  LCM  based 
gene  expression  profiling  yield  distinctly  different  lists  of 
metastasis-promoting  genes.  Contamination  of  the  tumor 
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cells,  and  cross  reactivity  of  mouse  RNA  to  human- specific 
chips  may  explain  these  differences,  and  suggests  that 
LCM-derived  data  may  be  more  accurate. 

Keywords  Breast  cancer  •  Estrogen  receptors  • 

Lymph  node  •  Lymphatic  vessel  •  Microarray  •  Metastasis  • 
Progesterone  receptors  •  Tumor  xenografts  •  ZsGreen 

Abbreviations 

CK18  Cytokeratin  18 

ER  Estrogen  receptor 

H  &  E  Hematoxylin  and  eosin 

LCM  Laser-capture  microdissection 
LN  Lymph  node 

Introduction 

Breast  cancer  metastasis  kills  over  400,000  worldwide  each 
year  [1].  At  diagnosis,  30-50%  of  primary  breast  cancers 
have  disseminated  to  other  locations  including  the  draining 
lymph  nodes  (LN).  Assessments  of  tumor  aggressiveness 
as  predictors  of  overall  survival  were  historically  confined 
to  histological  evaluations  by  trained  pathologists.  More 
recently,  molecular  approaches  have  been  shown  to  be 
superior  to  pathological  staging  to  determine  patient 
prognosis  and  treatment  options  [1-5]. 

Through  global  gene  expression  profiling,  breast  cancers 
have  been  shown  to  be  derived  from  at  least  two  func¬ 
tionally  distinct  cell  types,  basal  cells  and  luminal 
epithelial  cells;  each  responds  differently  to  therapeutics 
[6].  Lrom  these  genome- wide  analyses  new  hypotheses 
have  been  proposed,  which  suggest  that  the  ability  of  a 
breast  cancer  cell  to  metastasize  is  intrinsically  encoded 
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and  linked  to  the  same  genetic  alterations  that  initiated  the 
growth  of  the  primary  tumor  [7,  8].  Alternatively,  other 
data  indicate  that  only  certain  ‘rare  cells’  within  a  breast 
tumor  harbor  the  potential  for  tumor  growth  and  metastasis 
[9].  It  is  possible  that  both  views  of  progression  are  correct 
and  that  tumors  that  generate  metastatic  cells  differ  only 
slightly  in  their  transcriptome  from  the  bulk  of  the  cells 
that  comprise  the  primary  tumor. 

A  large  number  of  studies  have  been  conducted  to 
identify  the  genes  that  are  expressed  differently  in  primary 
tumors  compared  to  their  metastases,  and  that  may  there¬ 
fore  contribute  to  increased  aggressiveness.  Interestingly, 
the  published  ‘metastatic  gene  profiles’  are  not  composed 
of  the  same  gene  lists  [10].  Among  other  things,  these 
differences  among  studies  have  been  attributed  to  tissue 
sampling  variations  or  to  different  statistical  analyses  [10]. 

We  have  been  engaged  in  identifying  genes  that  may  be 
responsible  for  the  spread  and/or  growth  of  estrogen 
receptor  (ER)  positive  breast  cancers  to  the  LNs.  To  this 
end  we  determined  if  the  two  commonly  used  methods  to 
isolate  RNA — whole-organ  isolation  and  laser-capture 
microdissection  (LCM) — identify  similar  sets  of  genes  in 
an  experimental  model.  We  find  that  using  the  same  tumor 
lines  and  statistical  analyses,  gene  expression  profiles  from 
whole-organs  are  completely  different  from  the  profiles 
using  LCM  isolated  cells,  with  only  minor  overlap  between 
them.  Some  of  the  overlapping  genes  have  been  previously 
shown  to  be  involved  in  cancer  aggressiveness.  We  con¬ 
clude  that  the  methods  used  to  isolate  RNA  can  markedly 
skew  the  results  obtained. 


Materials  and  methods 

Fluorescent  xenograft  model  of  ER  +  LN  metastasis 

ER+  MCF-7  human  breast  cancer  cells  were  originally 
purchased  from  the  American  Type  Culture  Collection 
(Manassas,  VA).  The  generation  of  fluorescent  MCF-7 
cells  has  been  previously  described  [11].  In  brief,  ZsGreen 
retroviral  particles  (Clontech,  Mountain  View,  CA)  were 
isolated  from  PT-67  packaging  cells  (Clontech),  filtered, 
and  overlaid  onto  MCF-7  cells.  Cells  were  serially  trans¬ 
duced  and  then  subjected  to  G418  selection,  followed  by 
aseptic  fluorescence  associated  cell  sorting  to  isolate  a 
homogeneously  bright  green  subpopulation. 

All  animal  procedures  were  performed  under  a  protocol 
approved  by  the  University  of  Colorado  Institutional  Ani¬ 
mal  Care  and  Use  Committee.  Ovariectomized  athymic  nu/ 
nu  mice  were  obtained  from  Harlan  Sprague-Dawley 
(Indianapolis,  IN)  at  5-6  weeks  of  age.  To  establish 
tumors,  mice  were  anesthetized  with  Avertin  and  injected 
into  the  opening  of  the  lactiferous  duct  of  the  abdominal 
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mammary  glands  with  1  million  ZsGreen-expressing 
MCF-7  cells  in  100  microliters  of  100%  Matrigel  (Becton 
Dickinson,  Franklin  Lakes,  NJ).  Mice  were  also  implanted 
with  silastic  pellets  containing  17/f-estradiol  (2  mg  +  8  mg 
cellulose)  as  previously  described  [12].  Fluorescent  whole- 
body  imaging  (Illumatool  9900,  Lightools  Research,  En¬ 
cinitas,  CA)  was  done  weekly  to  determine  when  axillary 
LN  metastases  arose.  Mice  were  euthanized  by  C02 
asphyxiation  and  fluorescent  intravital  optical  imaging  was 
done  by  coupling  the  Illumatool  with  an  Olympus  (Mel¬ 
ville,  NY)  SZ-61  dissecting  microscope  and  Olympus  C- 
5050  digital  camera.  Fluorescent  tumors  in  mammary 
glands  and  matched  LN  metastases  were  removed. 

Whole-organ  RNA  isolation  and  expression  profiling 

Tumors  in  mammary  glands  and  LNs  containing  metasta¬ 
ses  were  extracted  and  frozen  at  -80°C.  RNA  was  prepared 
with  Trizol  according  to  the  manufacturer’s  specifications 
(Invitrogen,  Carlsbad,  CA).  Samples  were  in  vitro  tran¬ 
scribed,  biotinylated,  hybridized  to  Affymetrix  HG- 
U133+2  chips  that  interrogate  47,000  transcripts  with 
54,675  probesets  (Santa  Clara,  CA),  and  scanned  on  an 
Agilent  bioanalyzer  2100  (Santa  Clara,  CA). 

Laser-capture  microdissection,  RNA  isolation  and 
expression  profiling 

Mammary  gland  tumors  and  LNs  containing  metastases 
were  harvested,  placed  in  a  Tissue  Tek  cryomold  (EMS, 
Hatfield,  PA),  covered  with  Neg50  frozen  section  medium 
(Richard-Allan  Scientific,  Kalamazoo,  MI),  and  frozen  in 
isopentane  (Sigma,  St.  Louis,  MO)  cooled  by  liquid 
nitrogen.  Eight  micron  sections  were  cut  through  the  entire 
tumor  or  LN  with  a  cryotome  set  at  -20°C.  Each  section 
was  collected  on  an  uncharged  slide  (Fisher),  and  every 
10th  slide  was  stained  for  hematoxylin  and  eosin  (H  &  E) 
or  processed  for  cytokeratin  18  (CK18)  immunofluores¬ 
cence  to  confirm  location  of  the  cancer  cells.  Serial 
sections  of  slides  containing  cancer  cells  were  processed 
through  75%  ethanol,  water,  75,  95,  and  100%  ethanol  for 
30  s  each,  and  finally  dehydrated  with  xylene  for  1  min 
before  being  immediately  subjected  to  laser-capture 
microdissection  using  an  Arcturus  Autopix  (Sunnyvale, 
CA).  Two- thousand  cells  for  each  sample  were  extracted 
within  30  min  of  fixation  and  frozen  at  -70°C  until  further 
processing.  Suggested  protocols  from  Arcturus  were  fol¬ 
lowed  for  RNA  extraction  and  two  rounds  of  amplification 
using  the  Picopure  RNA  isolation  kit  and  RiboAmp  HS  kit 
(Arcturus),  respectively.  Samples  were  in  vitro  transcribed, 
biotinylated,  hybridized  to  Affymetrix  (Santa  Clara,  CA) 
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human  X3P  chips  that  interrogate  47,000  transcripts  with 
61,358  probesets,  and  scanned  on  an  Agilent  (Santa  Clara, 
CA)  bioanalyzer  2100. 

Statistical  analyses 

The  data  from  four  whole-organ  chips  (two  pairs  of  tumors 
and  their  matched  LNs),  and  four  LCM  chips  (two  pairs  of 
tumors  and  their  matched  LNs)  were  imported  into  Partek 
(Genomics  Suite  software  St.  Louis,  MO)  and  RMA  nor¬ 
malized.  For  comparative  analyses  between  the  whole- 
organ  and  LCM  data,  matching  probesets  from  the  two 
different  chip  types  were  identified.  Of  the  54,675  U133+2 
probesets,  3,487  were  dropped  from  the  analysis  because 
there  was  no  X3P  match,  but  7,301  probesets  were  used 
multiple  times.  The  result  was  that  all  but  62  of  the  61,358 
X3P  probesets  were  matched  to  a  U133+2  probeset.  The 
expression  values  for  each  organ  corresponding  to  either 
the  U133+2  or  X3P  probeset  were  exported  into  Microsoft 
Excel  (Redmond,  WA).  An  array-comparison  datasheet 
that  matched  the  U133+2  probesets  to  the  X3P  probesets 
was  obtained  from  Affymetrix.  Using  this  program  the 
RMA  normalized  expression  values  from  Partek  for  the 
four  U1 33+2  chips  were  matched  to  the  corresponding  X3P 
probesets  of  the  four  X3P  chips.  This  resulted  in  normal¬ 
ized  expression  levels  from  all  eight  chips.  The  chip  data 
were  also  imported  into  Genespring  (Agilent)  to  obtain 
present/marginal/absent  calls  for  each  chip.  Probesets  with 
absent  calls  on  all  eight  chips  of  the  same  type  were 
excluded  from  analysis  leaving  36,682  genes  that  were 
present  on  one  of  the  eight  chips  for  analysis.  These  genes 
were  imported  into  Partek  according  to  the  X3P  probeset 
identifiers.  Separate  paired  Mests  were  performed  on  the 
whole-organ  and  LCM  matched  primary  tumor  and  LN 
pairs.  P- values  <0.05  were  considered  significant. 

Gene  networking  analyses 

To  define  relationships  among  genes  that  were  differen¬ 
tially  regulated  in  tumors  versus  LN  metastases,  pathway 
analyses  were  performed  using  Ingenuity  software  (Red¬ 
wood  City,  CA).  Probeset  identification  numbers  from 
transcripts  with  increased  or  decreased  expression  in  LN 
metastases  were  used  for  gene  identification. 

Results 

Tissue  acquisition 

Estrogen-dependent  MCF-7  tumor  growth  is  required  for 
LN  metastasis  [11].  For  these  studies  MCF7  +  ZsGreen 


cells  were  grown  in  the  mammary  gland  of  ovariectomized 
immunocompromised  mice  that  had  been  supplemented 
with  an  estradiol  releasing  pellet.  Identification  of  axillary 
LN  metastases  was  achieved  with  fluorescent  whole  body 
imaging  (Fig.  la).  At  12  weeks  intact  mammary  gland 
tumors  (Fig.  lb)  and  their  matched  LN  metastases  (Fig.  lc) 
were  harvested  and  frozen  as  intact  whole-organs  or  cut 
into  8 -micron  frozen  sections.  LNs  and  tumors  were  then 
probed  with  a  breast  cancer  cell  marker,  cytokeratin  18,  to 
map  cancer  cell  locations  (Fig.  Id).  Adjacent  sections  were 
then  used  for  LCM  (Fig.  le-g).  The  LN  metastases  that 
were  used  for  the  intact  whole-organ  analysis  contained 
approximately  equal  amounts  of  cancer  cells  as  assessed  by 
ZsGreen  fluorescence,  CK18  immunohistochemistry,  and 
H&E  staining. 

Gene  expression  profiling 

To  compare  gene  expression  signatures  in  whole  tumors 
and  matched  intact  LN  metastases,  a  paired  Mest  was 
performed.  With  this  approach  1,930  genes  were  identified 
as  being  differentially  expressed  in  tumors  compared  to  LN 
metastases,  with  1,213  decreased  and  717  increased  in  the 
metastases  compared  to  the  tumors  from  which  they  orig¬ 
inated  (Fig.  2a).  For  the  LCM  samples,  the  same 
procedures  were  followed,  which  identified  368  decreased 
transcripts  and  913  increased  transcripts  in  the  metastases 
compared  to  the  tumors  from  which  they  originated 
(Fig.  2b).  Clearly  the  pattern  of  upregulation  and  down- 
regulation  differs  in  the  two  types  of  analyses. 
Additionally,  the  volcano  plots  show  much  larger  differ¬ 
ences  in  ‘fold-changes’  between  tumors  and  LN  metastases 
that  were  laser-captured  compared  to  the  whole-organ 
samples. 

Genes  that  change  in  whole-organ  compared  to  LCM 

When  whole-organ  and  LCM  paired  r-test  results  were 
combined,  3,181  transcripts  were  found  to  have  signifi¬ 
cantly  increased  or  decreased  expression  levels  in  tumors 
versus  LN  metastases.  Of  the  3,181  however,  only  68  (2%) 
were  identified  as  changing  by  both  analyses  (Fig.  3a),  and 
when  the  direction  of  change  was  calculated  for  these  68 
transcripts;  38  were  found  to  change  in  opposite  directions 
by  the  two  analyses  (Fig.  3a).  The  remaining  30  transcripts 
(~  1%  of  all  genes)  that  changed  in  the  same  direction  in 
both  the  whole-organ  and  LCM  analyses  are  listed  in 
Table  1. 

Results  from  the  paired  t- tests  were  compared  to  identify 
all  transcripts  that  were  differentially  regulated  in  either 
whole-organ  or  LCM  analyses.  Using  both  data  sets,  1,570 
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Fig.  1  Isolation  of  tumors  and 
LN  metastases  used  for  “whole- 
organ”  or  laser-capture 
microdissection  followed  by 
gene  expression  profiling,  (a) 
Fluorescent  whole  body  image 
of  an  immunocompromised 
mouse  with  bilateral 
MCF7  +  ZsGreen  tumors  and 
an  axillary  FN  metastasis,  (b) 
ex  vivo  fluorescent  image  of  a 
MCF7  +  ZsGreen  tumor,  (c) 

FN  metastasis,  (d)  cytokeratin 
18  immunohistochemistry  of  an 
axillary  FN  metastasis,  CK18 
red,  DAPI,  blue,  (e) 
hematoxylin  and  eosin  stain  of  a 
FN  identifying  the  metastatic 
MCF7  +  ZsGreen  cells  (*),  (f) 
combined  images  of  a  serial 
section  of  d  with  the  cancer 
cells  marked  for  FCM,  (g)  Post- 
FCM  image  of  captured  MCF7 
cells  from  e 


transcripts  were  found  to  be  downregulated  in  LN  metas¬ 
tases  (Fig.  3b).  Of  these,  1,202  were  unique  to  the  whole- 
organ  analysis,  357  were  unique  to  LCM  analysis,  and  11 
genes  exhibited  decreased  expression  with  either  method. 
A  similar  trend  was  observed  among  1,611  genes  that  were 
upregulated  in  LN  metastases,  with  698  being  whole-organ 
specific,  894  being  LCM  specific,  and  19  common  to  the 
two  analyses  (Fig.  3c). 

Ingenuity  pathway  analyses  were  performed  to  deter¬ 
mine  if  the  large  differences  between  the  two  types  of 
analyses  could  nevertheless  be  due  to  changes  among 
similar  functional  networks.  The  whole-organ  data  identi¬ 
fied  pathways  encoding  GM-CSF  signaling,  insulin 
receptor  signaling,  amyloid  processing,  Jak/Stat  signaling, 
and  sonic  hedgehog  signaling  that  were  the  most  different 
in  the  LN  metastases.  The  LCM  data  identified  completely 
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different  pathways  that  had  changed  the  most  in  LN 
metastases  including  PI3K/AKT  signaling,  PTEN  signal¬ 
ing,  synthesis  and  degradation  of  ketone  bodies,  oxidative 
phosphorylation,  and  keratin  sulfate  biosynthesis.  Clearly 
both  the  individual  genes  and  the  broader  pathways  that  are 
altered  in  tumors  versus  LN  metastases  are  different 
depending  on  the  assays  that  are  performed. 

Discussion 

Numerous  studies  from  independent  labs  have  shown  that 
there  are  gene  expression  changes  in  cancer  cells  within 
primary  tumors  compared  to  cells  that  seed  to  different 
metastatic  microenvironments  [8,  13].  These  studies  sug¬ 
gest  that  the  propensity  for  development  of  metastasis  is 
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Fig.  2  Identification  of  the  magnitude  of  gene  expression  changes  in 
tumors  and  LN  metastases.  (a)  Volcano  plot  showing  the  1,930  genes 
that  were  significantly  different  between  tumors  and  LN  metastases  in 
the  whole-organ  analysis,  (b)  volcano  plot  showing  the  1,281  genes 
that  were  significantly  different  between  tumors  and  LN  metastases 
by  the  LCM  analysis.  Lor  both  plots  the  magnitude  of  fold  change  is 
shown  on  the  x-axis  and  the  significance  is  shown  on  the  y-axis 

due  to  genes  that  are  altered  during  the  initial  tumorigenic 
process.  The  genetic  changes  in  metastases  have  been 
attributed  to  differences  required  for  successful  coloniza¬ 
tion  and  growth  within  vastly  different  metastatic 
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Fig.  3  Determining  the  overlap  of  the  genes  that  changed  in  the 
whole-organ  and  LCM  analyses.  Significantly  changing  genes  have  a 
decrease  or  increase  in  the  direction  of  the  fold  change  without 
specifying  a  cut-off  value,  (a)  Categorization  of  the  68  genes  that 
were  changed  in  both  whole-organ  and  LCM  analyses.  Shown  are 
genes  that  decreased  or  increased  in  LN  metastases  with  both  analyses 
or  those  whose  expression  changed  in  the  opposite  direction  in  each 
analysis,  (b)  genes  with  decreased  expression  in  LN  metastases,  (c) 
genes  with  increased  expression  in  LN  metastases 


microenvironments  [14].  Whether  the  cells  that  are  capable 
of  metastasis  are  ‘rare  cells’  or  ‘cancer  stem-cells’,  or 
reflect  any  cell  within  the  tumor,  remains  under  debate  [8, 
9].  Regardless,  determining  the  genes  that  are  differentially 
expressed  between  cancer  cells  growing  within  the  breast 
compared  to  those  growing  in  other  sites  is  crucial  to 
understanding  the  mechanisms  that  allow  cells  to  seed  and 
grow  at  distant  sites.  Identifying  these  genes  and  rendering 
them  therapeutic  targets,  is  an  essential  first  step  towards 
inhibiting  lethal  metastases. 
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Table  1  Genes  that  increased  or  decreased  in  LN  metastases  by  both  analyses 

X3P  PSID 

UniGene  ID 

Gene  symbol 

Whole-tissue  LC(LN/MGT) 

LCM  LC(LN/MGT) 

Genes  with  increased  expression 

in  LN  metastases 

in  both  analyses 

219529_3p_at 

Hs.64746 

CLIC3 

1.1114 

1.0972 

Hs2.429028.1.Sl_3p_at 

Hs.632184 

CARHSP1 

1.0548 

1.0498 

gll545730_3p_at 

Hs.l  12569 

GAN 

1.0492 

1.1879 

g4557696_3p_a_at 

Hs.99936 

KRT10 

1.0431 

1.2789 

208549_3p_x_at 

Hs.591909 

PTMA 

1.0376 

1.1448 

g45 075 7 4_3  p_a_at 

Hs.279594 

TNLRSL1A 

1.0373 

1.1685 

g388890_3p_a_at 

Hs. 89643 

TKT 

1.0285 

1.2193 

g9295191_3p_s_at 

Hs.584881 

C3orf28 

1.0242 

1.3023 

gl2804874_3p_a_at 

Hs.226390 

RRM2 

1.0222 

1.3370 

g4505038_3p_at 

Hs.1116 

LTBR 

1.0181 

1.1803 

Hs.l21017.0.Sl_3p_at 

Hs. 121017 

HIST1H2AE 

1.0171 

1.3161 

g79599 1 8_3p_a_at 

Hs. 381126 

RPS14 

1.0153 

1.2751 

Hs.2793 1 .0.S  l_3p_at 

Hs.592295 

GNPNAT1 

1.0083 

1.1926 

g7705770_3p_at 

Hs. 19385 

ABHD5 

1.0062 

1.3309 

Hs.27239.0.Sl_3p_at 

Hs.27239 

ZDHHC5 

1.0057 

1.1618 

201217_3p_x_at 

Hs.l  19598 

RPL3 

1.0050 

1.1517 

211073_3p_x_at 

Hs.l  19598 

RPL3 

1.0043 

1.1548 

g 1 3 623 7 08_3p_s_at 

Hs.l  19598 

RPL3 

1.0042 

1.1483 

g45 0197 6_3p_a_at 

Hs.512815 

AP3D1 

1.0028 

1.2395 

Genes  with  decreased  expression 

in  LN  metastases 

in  both  analyses 

Hs. 74034. 1  .S  l_3p_at 

Hs.74034 

CAV1 

-1.14927 

-1.21563 

Hs.248 12.0.S  l_3p_at 

Hs.472027 

CDS2 

-1.01957 

-1.06974 

Hs.  172329.  l.Sl_3p_at 

Hs.485628 

ZNL451 

-1.01761 

-1.06306 

Hs.l61962.0.Sl_3p_at 

Hs. 370510 

IGSL4 

-1.03359 

-1.01245 

Hs.  145612. 0.S2_3p_at 

Hs.  190622 

DDX58 

-1.05818 

-1.26471 

Hs.  106440.0.  Al_3p_at 

- 

- 

-1.07679 

-1.05797 

g9951924_3p_at 

Hs.210995 

CA12 

-1.08361 

-1.05245 

g7661679_3p_a_at 

Hs.455336 

NELL 

-1.01845 

-1.04308 

g4929370_3p_x_at 

Hs.465885 

ILL3 

-1.0444 

-1.0716 

gl2052731 _3p_x_at 

Hs. 35 1906 

ZNL506 

-1.04132 

-1.0935 

gl 205273 l_3p_at 

Hs. 35 1906 

ZNL506 

-1.04186 

-1.07409 

Shown  are  the  probeset  identification  code  (PSID),  the  UniGene  identification  code,  the  gene  symbol  and  the  RMA-normalized  fold  changes 
from  the  whole-organ  and  LCM  analyses 


For  the  whole-organ  analysis,  RNA  was  isolated  from 
entire  tumors,  or  from  whole  LNs  containing  similar 
amounts  of  cancer  cells  within  them.  For  LCM,  RNA  was 
isolated  from  2,000  pure  cancer  cells  of  tumors  and  LN 
metastases  (Fig.  1).  To  compare  these  sets  we  chose  an 
approach  that  is  commonly  used  to  study  large  numbers  of 
gene  expression  changes;  namely,  global  gene  expression 
profiling  with  Affymetrix  gene  chips.  The  chips  used  are 
designed  to  theoretically  hybridize  only  with  human  RNA, 
and  therefore  both  the  whole-organ  and  LCM  approaches 
should  inform  us  only  of  the  gene  expression  changes  that 
occur  in  the  human  cancer  cells  within  the  mouse  host 
sites.  It  was  therefore  surprising  to  find  greater  Told 
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changes’  in  LN-specific  genes  assessed  by  LCM  (Fig.  2a 
vs.  Fig.  2b). 

Since  the  gene  chips  used  for  the  two  methods  interro¬ 
gated  the  same  transcripts,  with  all  but  62  of  the  61,358 
X3P  probesets  having  a  matching  U133+2  chip  probeset,  it 
was  surprising  that  very  few  genes  changed  similarly 
between  the  two  analyses  (Fig.  3).  Furthermore,  of  the  68 
transcripts  that  were  significantly  different  in  LNs  by  both 
analyses,  38  were  regulated  in  the  opposite  direction,  being 
increased  in  the  metastases  by  one  analysis,  and  decreased 
by  the  other  (Fig.  3a).  These  results  suggest  that  the  Af¬ 
fymetrix  chips  are  binding  not  only  human  RNA  but  also 
mouse  RNA.  This  aberrant  hybridization  may  be 
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responsible  not  only  for  the  gene  differences  reported  by 
the  two  methods  but  also  for  the  lack  of  overlap  between 
the  functional  pathways  identified. 

No  other  reports  have  determined  how  sample  collection 
methods  influence  LN  metastasis  gene  signatures.  How¬ 
ever,  one  study  [15]  compared  28  breast  tumors,  17  ER 
positive  versus  11  ER  negative,  and  identified  genes  dif¬ 
ferently  expressed  between  these  two  cancer  types  using 
both  LCM  and  bulk- tumor  gene  expression  profiling. 
Interestingly,  of  all  the  genes  found  to  be  significantly 
different  between  ER  positive  and  ER  negative  tumors,  less 
than  one-third  were  common  to  both  analyses.  These 
findings  are  interesting  because  the  same  human  tumors 
were  used  for  comparison,  yet  the  ER  signature  changes 
were  dependent  on  the  type  of  analysis  conducted. 

We  suggest  two  possible  explanations  for  the  differ¬ 
ences  between  our  whole  organ  and  LCM-generated  data. 
(1)  It  is  possible  that  multiple  unique  transcripts  are  altered 
by  the  amplification  procedure  required  for  LCM.  How¬ 
ever,  a  previous  study  showed  that  amplification  from  the 
3'  end  of  transcripts  does  not  have  a  major  impact  on 
overall  gene  profiles  with  Affymetrix  GeneChips  because 
the  probe  sets  on  the  arrays  are  designed  using  3'  end 
sequences  [16].  Additionally,  another  report  found  no 
significant  differences  in  the  expression  intensity  of  21 
housekeeping  genes  from  matched  LCM  and  bulk-tumor 
gene  expression  data  [15].  (2)  An  alternative  hypothesis 
posits  that  Affymetrix  chips  can  bind  not  only  human 
RNA,  but  also  mouse  RNA,  thus  altering  the  chip  readout 
from  whole  organ  analyses.  Naef  et  al.  [17]  recently  tested 
this  idea  and  found  that  indeed,  mouse  RNA  hybridized  to 
the  “human- specific”  gene  chips.  Interestingly,  we  found 
that  on  average  16.4%  more  probesets  were  present  on 
whole-organ  samples  compared  to  LCM  samples,  sug¬ 
gesting  that  in  xenografts,  mouse  RNA  was  contaminating 
the  whole-organ  analysis.  The  studies  by  Naef  et  al.  sug¬ 
gest  not  only  the  possibility  of  cross-species  contamination, 
but  the  possible  skewing  of  data  due  to  intraspecies  con¬ 
tamination  of  tumor  cells  by  normal  epithelial,  stromal, 
vascular,  and  immune  tissue  elements  within  an  organ. 
Thus  it  may  be  critical  to  use  LCM  to  determine  how  each 
tissue  component  of  an  organ  contributes  to  the  overall 
gene  signature.  Based  on  this,  we  suspect  that  our  LCM 
data  represent  a  more  accurate  gene  expression  comparison 
between  tumors  and  their  matched  LN  metastases. 

These  studies  nevertheless  identified  30  transcripts  that 
were  increased  or  decreased  in  LN  metastases  by  both 
analyses  (Table  1).  These  genes  may  represent  the  best 
targets  for  further  studies  to  identify  the  mechanisms  that 
contribute  to  the  spread  or  growth  of  breast  cancer  cells 
into  and  within  LNs.  Prothymosin  alpha  (PTMA)  was 
increased  in  the  LN  metastases  and  has  previously  been 
shown  to  increase  proliferation  rates  of  breast  cancer  cells 


[18].  This  finding  is  interesting  in  view  of  our  previous 
findings  of  increased  proliferation  rates  in  LN  metastases 
[11].  Tumor  necrosis  factor  receptor  superfamily  member 
1A  (TNLRSL1A)  was  also  increased  in  the  LN  metasta¬ 
ses.  Its  expression  has  previously  been  shown  to  decrease 
apoptosis  in  MCL-7  cells  treated  with  human  TNL  alpha 
protein  [19].  Ribonucleotide  reductase  M2  polypeptide 
(RRM2)  and  lymphotoxin  beta  receptor  (LTBR)  were 
both  increased  in  the  LN  metastases,  and  both  increase 
cancer  cell  survival  [20,  21].  Interestingly,  RRM2  is  tar¬ 
geted  by  gemcitabine  and  hydroxyurea  in  lung  and  skin 
cancers,  respectively  [22].  A  transcript  whose  expression 
was  decreased  in  LN  metastases  and  has  been  studied  for 
its  effects  on  cancer  progression  is  caveolin-1  (CAV1). 
This  scaffolding  protein  is  the  main  component  of  the 
caveolae  in  plasma  membranes  of  most  cell  types,  and  its 
expression  inhibits  proliferation  in  human  breast  cancer 
cells  [23].  Its  decreased  expression  may  therefore  play  a 
role  in  the  increased  cell  proliferation  seen  in  LN 
metastases  [11,  24]. 

In  summary,  this  report  evaluated  whole-organ  versus 
LCM  for  gene  expression  profiling,  in  order  to  identify 
therapeutic  targets  that  may  reduce  breast  cancer  spread  to 
and  growth  in  LNs.  However,  when  comparing  these  two 
experimental  approaches,  we  found  almost  no  overlap 
among  the  transcripts  and  functional  pathways  identified. 
The  30  transcripts  that  were  common  to  both  analyses  may 
be  the  most  promising  for  future  therapeutic  targeting. 
However,  our  data  clearly  demonstrate  that  the  methods 
used  to  prepare  tumors  for  gene  expression  profiling  can 
markedly  skew  the  results  obtained. 
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